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1 
» RAPPELS SUR LA CINEMATIQUE DES FLUIDES ET LA 


DYNAMIQUE DES FLUIDES INCOMPRESSIBLES 


>» 1. Equation de continuité (conservation de la masse) 
| | L'équation de continuité doit traduire le principe de conservation de la masse т: 


la variation de masse pendant un temps dt d'un élément de volume fluide doit être 
+ égale à la somme des masses de fluide entrant diminuée de celle de fluide 
sortant. 


Bilan de masse sur dv = dx,dx,dx; pendant l'unité de temps dt : 


Xs u,(x;+dx,) 
(о u, dx,dx;dt |. — Lo u, dx; dx. di], -— )= (о и, l. — Lo u, l. idi )dx,dx,dt 
u,(X;+dx;) lo Ua dx, dx; di |, i lo и dx, dx; dl +dx, ) m [> Hj L. n [р H: Ke Jax, dx; dt 


(о u,dx,dx,dt]. -|p u;dx,dx;dt|, 4, j= (о nl Jon, )ax, dx, dt 


из(хз) 
u,(X, dx) 
Développement en série de Taylor à l'ordre 1 de lo и; | „ax puis calcul de la différence: 


lo u; 1 +dx, =[р и; l. m M Au la, Ld L Hi 1 -l Hi l. +ах Ја, dt =— дри) ln dx,dx;dt 


i 
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` Au total, le bilan à travers les 6 faces donne : 


4 


h x, OX, í. 


S|- Agu) dx, dx,dx,dt => Е t lor. ` == Tam 
X; 
> Or la masse élémentaire dm est donnée par : 


dm — А ах,ах,ах, dt — — = CS 1 E ах,ах,ах, gf 
Kë ^ а X5 TA 


C'est l'équation de la continuité, qui traduit la conservation de la masse. 


| ou encore 


Z 3 
| 
RAPPELS SUR LA CINEMATIQUE DES FLUIDES ET LA 


` ` DYNAMIQUE DES FLUIDES INCOMPRESSIBLES 


1 Toute dérivée par rapport au temps est nulle: d... =0 


i Ot 


і AP) ai (рӯ)- div(pV )= p.div(V }+V.grad(p)=0 


SW Ecoulement isovolume (écoulement incompressible) 


` p-cste "mj raid (P) = 0 —— div(V )= 0 
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7 2. Accélération de l'écoulement 
» Selon la définition de la dérivée particulaire pour un vecteur, au sens cinématique, le 
vecteur accélération d'une particule fluide dans l'écoulement s'écrit : 


…._dV _ôV ду Ôx; OV ду âx, OV Ox, OV ôx; 
са Or ôx; Or а ôx, Ot, Ox, Ot, Ox, Gt, 


=u; P =и; 
„© F ii Ou, + Ou, 
И р 1 ôx, Ох» 0х; 
OV OV Ou» Ои» ди» 
U3 — = > + u] + И? — `+ U3 === 
OX; Ot OX] Ох» 0X3 
диз 2 диз 3-3 диз 
1 TE dE 
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> Ou encore : Ou, Ou, Ou, 


oV ди, ди, ôu, и, ZE E bk Ji 


p=% 
Ot | Ox, Ox, Ф, Ot 
Ox, Ox, Or 


… Le tenseur gradient est donné par : 


ex, OX, ox, 


grad )-v(e)- Z de M (A) u) 


Á 
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> RAPPELS SUR LA CINEMATIQUE DES FLUIDES ET LA 


DYNAMIQUE DES FLUIDES INCOMPRESSIBLES 


> 3. Forces surfaciques - Tenseur des contraintes 


». Pour un fluide réel (visqueux) en mouvement, les forces de surface ne sont plus 


seulement normales à la surface : il existe également des contraintes tangentielles 
. qui sont dues à la viscosité du fluide (frottements). 


. En un point M d'une surface élémentaire dS, la force de surface est donnée par : 


| 


l 


dF =T. dS 


avec Т est la contrainte s'exerçant sur la 
surface élémentaire de normale 77. 
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> Considérons une surface | à l'axe x,. La normale à cette 
. surface est : п = e. 
1 


La contrainte exercée sur cette surface, notée D , peut 
1 
~ être décomposé ainsi : 


— 


27] 


> Les composantes с, et g,,sont des composantes 
. tangentielles : on les notera par ç, et 7,, pour les 


distinguer de la composante normale с, . Xa 
On... 


_ Considérons la sur face | à l'axe х,. La contrainte 

s'écrit : " 

T = È б z €» 
= T,2€, T O55€; + 15€; 


X2 


_ Pour la sur face | à Гахе x,. La contrainte s'exprime : © 


T. —T,4€, +T,;@, +033€3 


X3 
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— Pour une surface d'orientation quelconque, la normale peut se décomposer dans 
. le repère cartésien ainsi : 


n=n,e, +п,е, tne, 

. La contrainte s'exercant sur cette surface s'écrit : 
L =n;T. +n, +n. 

… Aprés développement, on obtient : 


el, =n, o, ë, +7, + Te; )+ n; (2,2, +0,,е, +T32C3 )+ П; (ЖЗ +т›уе,+о;;е; ) 


— 


—VtOj; tnt nz, + (n, Ty t 055 DAD + (n, 7311,95 t N3033 Je, 


» En écriture matricielle, cela revient à : 
От T; з || ñi 


1,=| fan Gas Tz ||һ, 


T3; T32 Озз )\ 113 Tenseur des contraintes 
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Le tenseur T peut étre décomposé en deux parties : 

E c, 0 

- Un tenseur sphérique umm» 0 c. 
0 0 

- Un tenseur de trace nulle [== "7 т, 


1 Oz 


forcenormale 
à la surface 


L Force de pression hydrostatique : œ = -p 
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» 4, Principe Fondamental de la Dynamique 
» Nous allons établir l'équation fondamentale de la dynamique en appliquant le 
principe fondamental de la dynamique à un élément de volume fluide en 
Ku n. Nous ferons par la suite un bilan des forces s'exerçant en surface et 


B. volume. 


dF = UE, + dE, = pd v — 


Forces de Forces de (e N 0 


dF, ë, edF, ë, +dF; ё, 


L Prenons par exemple cette composante suivant l'axe x.. 
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Jar, sel, kl äs 
T (с 22 e = lc 22 l Јах,ах, PRO 


ш (г. Ka ü D l. )dx,dx, 


Développement en série de Taylor à l'ordre 1 : 
Or 

D = D) +— dx, 
Ox, 


im) kl, E lol, + LP dx, 


OX, 


DAMM = rl. + — dx; 
X3 
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DYNAMIQUE DES FLUIDES INCOMPRESSIBLES 


Finalement, on obtient sur l'axe x; : 


De la méme facon, on obtient sur les autres axes x, et x, : 


шщ) AXCS X, 


umm AXCS X; 


OT. 

Se Ou Tp 

D'où l'expression générale "ME = = = а= т о, 
X 

дс. Zar T32 


Үү HY HY A = 
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Or, l'effort élémentaire total dF est la somme globale des efforts dus aux 
contraintes surfaciques dFs et des efforts dus aux contraintes volumiques dF, : 


dF = dF, + dF, na 


4 v 


Dans le tenseur des contraintes, en faisant apparaitre les deux parties symétrique 
et antisymétrique, on obtient : 


T = pl + T'u ç VT = 4 


> о ҹу 


0 0 8 = 
p 0 +VT' --Vp4VT' 
0 p 


C'est l'équation fondamentale de 
la dynamique, écrite sous une 
forme locale 


Finalement : 
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E 
a шш 
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ч 


> Equations de Navier-Stokes pour un fluide newtonien 


». Fluides newtoniens > les composantes du tenseur des contraintes visqueuses 
T' dépendent linéairement des composantes du tenseur des taux de déformation 


Í pure e. 
> Normalement, une rotation pure n'engendre aucune déformation паа pas de 
contrainte: c'est-à-dire 7' еї œ ne sont pas liés. 


E | | E E | ди | 
à. éléments du tenseur <= sont: £= m EE + 2) => Тепѕеиг symétrique 
X. X. 
J I 


2 


Pour un fluide isotrope, les éléments де 7' еї с sont liée par l'expression : 


Oj =2 HU Ey + // (Ж T€ +; )ó; 


> Ou, = 
avec: Z, +E, +E33 = LS y, = div [V ) 
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Pour un fluide incompressible newtonien, l'équation fondamentale de la 
„dynamique devient : 


РИМЕ fa cei E м 


dt 5 
>. = у= &.. . Ou. 
où le terme Vzs'écrit : V ç = Kä >` P dc wasapa E 0 
| d| е. J ° * Ql; Ox 
Ce qui donne : WE j Í ' 
D - 1( d'u. Ou. ] Он. LEE: ди. 
VE m — - + É е. = — i e; T— uq - e; 
— Ox; ) a B = US > 3 
| I "ee — l car le fluide est 
Hi ` incompressible 
p V(V.v) 
D'où l'équation fondamentale de la dynamique VES AV 
pour un fluide newtonien incompressible : 2 


i imp C'est l'équation de Navier-Stokes 
(équation locale) 
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Dérivée particulaire de la vitesse y: 


dV EC OV qx) on ol 
dt Ot 
Dérivée srivée 
temporelle convective 
En posant g — — 8 е, l'équation de Navier-Stokes projetée sur les 3 axes du 


repère cartésien, (о,х,,х,,хз), s'écrit alors : 


Ô UU, C u Ou, 
H 2 = 2) E E 
ОХ, ex OX; 


= Ә?и, а | 
+ Zí 2 ЕХ 2 nb 2 
ОХ; Cx OX; 


= Ә?и, | 
+ и 3 + + — og 


2 2) E 
ОХ, ОХ, ОХ 


Ces équations déduites du principe fondamental de la dynamique pour un fluide newtonien 
incompressible, traduisent la conservation de la quantité de mouvement. 


АШ HY HY A = 


м 
18 
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Cas d'un fluide parfait incompressible 


Un fluide parfait est un fluide au sein duquel on peut négliger les effet visqueux, 
c'est-à-dire les contraintes liées à la viscosité. Ainsi, pour un écoulement 
incompressible les équations de Navier stokes deviennent : 


а 27 rel en 


dt (а 


W HY HY A LA 


19 
RAPPELS SUR LA CINEMATIQUE DES FLUIDES ET LA 


DYNAMIQUE DES FLUIDES INCOMPRESSIBLES 


L'équation de Bernoulli traduit la conservation de l'énergie de fluide en 
écoulement sur une ligne de courant, dans un écoulement permanent de fluide 
parfait méme en présence de tourbillons. Ainsi, la somme des énergies cinétique, 
potentielle et celle communiquée par les forces de pression se conserve. 


C'est-à-dire, entre deux sections droites de l'écoulement monodimensionnel, 


notées par exemple 1 et 2, on obtient : 


Il va de soi que si, entre les deux sections 1 et 2, l'écoulement recoit ou perd une 
énergie W par unité de masse, on aura alors : 


Cette équation s'appelle équation de Bernoulli généralisée 
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Cette équation de mouvement pour un fluide parfait incompressible, s'écrit 
également : 


А _ В у 
p dV --Vptpg TE TA = BEEN 


— V? mM _ _ 
IT, = + ZE) + rot(V )^ V = Së grad(p)+ g 
p 


E. 


Er 


+ rot\V )A V = к grad(p)+ g- 2L wins grad(- g х, ) 
p p 
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| 


Pour un écoulement monodimensionnel, rotationnel ou irrotationel, les équations 
d'Euler le long d'une ligne de courant, aprés projection sur 4/ se réduisent а: 
у? E > 


srad| as = — ^ (zrad( p))d? — (rad (e x; lui 


| car 


— 


rot(V )A V est perpendiculaire à V, et V est colinéaire à dé 


Puisqu'il n'y a qu'une coordonnée de l'espace, l'équation précédente devient si on 


v V? 
CS are L P e+ g LD ae d |r Leg =0 
p p 


"Aprés intégration de cette équation différentielle du 1* ordre à variables 


. Séparables, on obtient : 
p*ipV spezcOn formule de Bernoulli (1738) 


4 Vi р 
1 Sa e + 27 = Cste OU encore : 
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— 5. Théoréme d'Euler ou théoréme des quantités de mouvement 


Le théoréme de Bernoulli est d'une utilisation trop limitée. En effet, il ne permet 

pas d'exprimer les actions mécaniques pouvant apparaître entre des fluides et 
_ des solides par exemple. D'où la nécessité d'introduire un deuxième théorème: 
c’est le théorème d'Euler. 

Ce théorème est établi à partir de la relation fondamentale de la dynamique pour 

un fluide parfait (incompressible, non visqueux) et un écoulement stationnaire : 


dérivée particulaire 


forces de pression forces volumiques (poids) 
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Application à un tube de courant : 

Appliquons le cas précédent à un tube de courant dans lequel la vitesse dans une 
section droite est uniforme. L'écoulement d'un fluide parfait esas, est 
considéré comme étant permanent. 


R + P = [а= (ла + ffias + (7) та 
mer li, fares 5328.48.48, S, | S, | S U 


- pV, V, S, car V=Cste sur S, pv, V. S, car V=Cste sur S, =0 car V Lñ 


Chapitre 1 


Ecoulements compressibles d'un gaz parfait et idéal 


4 Т 


Chapitre 3 : ET 


DYNAMIQUE DES FLUIDES COMPRESSIBLES 
ÉCOULEMENT MONODIMENSIONNEL EN RÉGIME 
PERMANENT D'UN FLUIDE COMPRESSIBLE (GAZ PARFAIT) 


1. Introduction 


L'hypothèse de fluide incompressible est évidemment parfaitement justifiée pour 
les liquides. Dans le cas des gaz, elle est encore valable tant que les variations de 


œ pression mises en jeu restent modérées mais toujours petites par rapport à la 


pression absolue du fluide : les calculs s'effectuent alors en supposant le milieu 
isotherme (c'est ce que nous avons fait jusqu'à présent). 


Mais d'importantes variations de pression entrainent des variations notables de la 
masse volumique p, et on ne peut plus considérer dans les calculs o comme une 


—— constante. 


De telles variations importantes de pression interviennent notamment dans le cas 


œ~ des grandes vitesses (approchant la vitesse du son) d'écoulement d'un gaz 


(tuyère) ou dans celui des déplacements rapides d'un corps solide dans un milieu 
gazeux (avion, lancement d'un projectile, turbine à vapeur). 
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ÉCOULEMENT MONODIMENSIONNEL EN RÉGIME 
PERMANENT D'UN FLUIDE COMPRESSIBLE (GAZ PARFAIT) 


> 2. Relation de Saint- Venant 
Pour établir la relation de Saint-venant, nous allons considéré les hypothèses 
suivantes pour le gaz : 


- Écoulement permanent 

- Écoulement irrotationnel 

- Écoulement monodimentionnel selon une abscisse curviligne l 

- Écoulement soumis à la seule force à distance qui est celle de la pesanteur 

- Fluide supposé parfait lors de l'écoulement (sans viscosité, ou sans frottement ou sans énergie dissipative) 
- Écoulement adiabatique 


Ces hypothèses conduisent à l'établissement de la relation de Saint-Venant, qui 
. est bien vérifiée expérimentalement pour l'air en écoulement transsonique ou 
supersonique. 
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| Compte tenu de ces hypothèses, l'équation de Navier - Stokes se réduit à 


> l'équation а Euler : s^ 
NM 
EE es)v)-otsye-»(Iste)-v^fsv))-- So p 
ou bien 


(| V) er s " E 
о) 6л) S oa So o x) 
~ ou bien: 
V progx,+Zpv?)=pV^(7 ^Y) 
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> Pour un écoulement monodimensionnel, rotationnel ou irrotationel, les équations 
d'Euler le long d'une ligne de courant, après projection sur 4/ se réduisent à : 


—U 


ed |й = = (erad(p)}di —(grad(s x, ))a? " eo) db) 
rotlV |A ui 


— 


car 
rot(V )^ V est perpendicu laire à V, et V est colinéaire à а? 


C 

* 
e 

ei 

ei 

EE EE Ke NA VV TT Lis 
BM БАБАЕ 
ei 
À 
wë 


Puisqu'il n'y a qu'une coordonnée de l'espace, l'équation précédente devient si on 


« pose X, = Z ` 

y? ү? 
SE — dta L P qe 4 Ou-a — u C nc 
dí! 2 p dt dt 2 p 


Dans cette équation différentielle du 1% ordre, la masse volumique р n'est pas 
constante. Donc, cette équation ne peut étre intégrée directement. 
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> D'où: 


2 


V 
d| — EET 
2) p 


. Dans cette équation, on apercoit les termes d'énergie cinétique et d'énergie 


| potentielle du gaz qui découlent directement du premier principe généralisé de la 
thermodynamique établi pour un système ouvert. 


— L'enthalpie H, est définie par : 


H —U-- pw 

. Sa variation élémentaire dH s'écrit : 

| dH = 40  pv)- dU + p dv +v dp = óW + 80 + p dv +v dp 
——p dv +00 + p dv +V dp = дО +V dp 


— Donc: 


dH =óQ+v dp = óQ4 dp ou bien de vpn E 
p p 
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. Soit: 


2 


Vg V —— 
5 P gdz = Ж — 80 )+ ва: =0 
2] p 2 
. Comme l'écoulement de fluide est supposé adiabatique, donc 50 — 0. D'où : 
V E 
[5 em + gdz = 0 


> Après intégration, on obtient : 


Cette équation traduit la conservation de l'énergie d'un fluide en écoulement 
> (énergie mécanique et énergie thermique (enthalpie massique)), quand on ne lui 
apporte aucune énergie extérieure : ni énergie mécanique, ni chaleur. 
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Comme par hypothèse, le fluide est supposé parfait lors de l'écoulement (sans viscosité, 
ou sans frottement ou bien sans énergie dissipative), alors dans ce cas опа: 


H -C,T + Cte = 4 rT + Cte =" pV + Cte 2 —— P. + Cte 
y -1 pest jp 


. D'où l'équation précédente devient : 


"mmm Équation de Saint-Venant 


` Cette relation est à la base de toute l'aérodynamique des écoulements à grande 
vitesse. En général, on peut négliger le terme d'énergie potentielle du à la 


. pesanteur, et on écrit : 
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La relation de Saint-Venant peut étre également obtenue à partir de l'expression 
d du premier principe généralisé de la thermodynamique pour un systeme ouvert. 
| En effet, pour un écoulement en régime permanent, on a : 


2 2 
=) 2, Hut oar -m, EETA -0 
dt Volume de Contróle dt dt 2 2 


Volume de contróle 


€ д — 
Н, di 


Dans cette relation, P, et représentent respectivement les puissances et 
mécanique et calorifique apportées au volume de contróle pendant l'unité de 
temps, m, et m, sont les débits massiques entrant et sortant de ce volume. 


+ Si le régime est permanent mp, 7n, = m. 
* Si l'écoulement est adiabatique et oW 
| | k |u — = 20 = 0 
s’il n’y a pas de travail de frottement dt dt 
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ÉCOULEMENT MONODIMENSIONNEL EN RÉGIME 
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> Ainsi, l'expression précédente obtenue à partir du premier principe de la 
thermodynamique généralisé pour les systémes ouverts s'écrit : 


2 2 


. Cette relation est trés générale. Comme nous l'avons vu précédemment, c'est elle 
qui conduit, en fluide incompressible, à la formule de Bernoulli (H — p! p, + Cste ). 


_ + Si l'écoulement du gaz parfait est adiabatique et de plus réversible (ce qui est 
. l'hypothèse généralement utilisée), les lois d'adiabaticité réversible : 


pV” = Ё. = Cste 
p 
Iu 


p'T — Cste 
Sont alors aussi applicables. 
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» 3. Application aux écoulements supersoniques 
3.1. Notions sur les écoulements supersoniques 
Les écoulements supersoniques sont des écoulements à grande vitesse d'un gaz 
`>». ou d'une vapeur, fluides pour lesquels les effets de la pesanteur peuvent le plus 
souvent étre négligés. 
Dans le cas de l'air et des gaz de combustion, l'expérience montre que l'hypothèse 
~ d'écoulement adiabatique réversible d'un gaz parfait est bien vérifiée. On se 
limitera alors dans notre cours à l'étude de ce cas qui correspond aux souffleries 
supersoniques et aux lancements de fusées. 
`> Les relations liant la pression p, la masse volumique p, la température T et la 
vitesse V sont alors : 
+ L'équation d'état du gaz parfait mmm p/p-(R/M)T —rT 
+ La relation d'adiabacité réversible (isentropique) 1} P. — Cste 
+ La EE de A Venant écrite sous ses deux formes équivalentes : 


E — |Z = Cste OU bien 
2 Y — JI P 


Vv С 
+C T = Cste 
2 P 
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~ Comme la vitesse V de l'écoulement de gaz est voisine de celle du son, il est alors 
important de la déterminer en se référant à la propagation d'une onde plane dans 
une colonne de gaz. 


x x + &(х) x + dx x + dx + e dx) 


à l'instant t à l'instant t+ dt 


Le principe fondamental de la dynamique s'écrit : 


" / e 
> F =m £= (es dx) E = ((р, +æ} (p + ах) 5 = 2200 суз 


dé alp) 
—Ó93 2. > 
La propagation d'une onde plane s'effectue sans échange de chaleur, c'est-à-dire 


de manière adiabatique (réversible). La relation entre la variation de pression óp 
(dp) et la variation du volume (dv) de la tranche de fluide est déduite de pv = Cste 
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P SE 
dv dv p d d 
n> óp = dp=-y р --y р = peter p o| 22) р 


or. v=Sx > dv =S dx =S ((x+dx+ а) (+ (0) на) 8 € 
X 
. dv 0€ дё SE ЖА 
d : = = et = dp = -y p— E yyy 
ñi Md V Ох š m a PX d Ох ‚аша 
CN: 


2 
- $9) = "n: Z bien = 
E on дг? 


д? 
Comme опа : 5 — = 
~ d T p ox 


>» On pose : a= |/ P шшщ); vitesse de propagation de l'onde (vitesse du son) 
О д? д? 
Donc : =а? 


ч : =a 
| д? дх? 
р dp dp dp) р 
отте : ру =Cste > = =Cste > —-y—-0 > bun ерй ü "e 
A p p p dp); p 


isentropique 
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`> On appelle nombre de Mach de l'écoulement, le rapport : 


vitesse locale de l'écoulement de gaz 
vitesse locale du son 


V Vv > 
M = = lee 1 MES 
SiV«a (M « 1) = um ue. Écoulement subsonique 
SiV=a(M=1 mm = Écoulement sonique 
Si V>a (М> 7) u = c= Écoulement supersonique 


`< Selon l'équation de Saint-Venant, la température d'un écoulement isentropique 


évolue avec sa vitesse. Ill en est de méme pour la vitesse du son a: 
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— Existence d'une vitesse limite pour le gaz (supposé parfait et en écoulement isentropique) : 


p, p, T, V 
| Pis Pi, T;, et V,= 0 sont les conditions génératrices 


Fluide parfait en écoulement isentropique шшш) pV” = = —Cste ou р” T —Cste 
y ү? P y? 
> Saint-Venant ess == di LX e и + Y Pi V Pi OU — + CT = С.Т, 
| 2 y-1p 2 y-1p; yeu 9 
Si a est la vitesse du son, on écrit: q= ЁР _ JyrT 
p 
< y? 2 
Donc, Saint-Venant devient : V? + 2 y — 2 а; Il —+ 2 127. s: 
y —1 y—1 a у-1 y-la 
ч 2 Z T 
. Or, le nombre de Mach Mest : M =k шш)» M +⁄ = = < і 
a Y=] eat T 
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b =} I p ЕЯ 
Fluide parfait en écoulement isentropique imp p ” T = p,” T, ou = 


.. Or, le théorème de Saint-Venant s'écrit également : 


V =_= C — 
стл, => V?=2C,(T. -Т)=2 (r -T)-2 ry T. E 


y-1 
V? =g ZÉ д SEA t 
ye D, ye 


Sh 


q C'est la vitesse Vétablie en fonction des conditions génératrices de l'écoulement 
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> Cela veut dire que pour un écoulement où la détente pourrait être poussée jusqu'à 
une pression nulle (p = 0), la vitesse limite V, en ce point serait : 


.. Pour le cas de l'air, par exemple, nous avons : 


R _ 832 
M. 2910? 


air 


y-140et r = ~ 286,89 J kg KI 


— Ainsi, pour T = 20 °C = 293 K, ona: 


pu еу у LI Uus 
y-1 140-1 
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. 3.2. Écoulement monodimensionnel permanent dans une tuyére 
3.2.1. Théoréme de HUGONIOT 
Quelle relation liant la vitesse de l'écoulement de gaz dans la tuyëre à sa section ? 
… Soit l'écoulement permanent dans le canal ci-dessous de section S variable. 


S 


. Le débit massique à travers la section S est : m = p S V 


_ Comme ce débit massique se conserve, alors en différentiant son expression, on а: 
| dm d as av 
—— = KAN quee 
Or, en l'absence de forces à distance, on a : 
: V? dP dP 
a]l p- i 


d 
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Mais comme l'écoulement est isentropique, on peut remplacer dp par son 
expression en fonction de do obtenue de la valeur de la vitesse du son a 


2) 5 E _ 2 

a = | 22.) -—— аг -| © Donc: dp = а d 
e s dp js d P 
Ainsi, l'équation établie précédemment devient alors : 


vavi P уау EP = ou bien dp 
P p 


Or. = dp ds, dV ow ; -f dV + ds ,4V _ 
m p S V s 


2 
ou encore: d$ , dV GIE 


| a Soit : ГЕНИ, "a C'est la formule de HUGONIOT 
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~ Formule de HUGONIOT 


dS/dV est négatif : /a vitesse 


i V < < 1 (écoulement subsonique F=» | A ка 
3 a ou M<1( ade) » 1 si la section diminue 


(dS/dV) = 0 : la section est 
Si V =a ou M = 1 (écoulement sonique) n 4 minimale lorsque la vitesse V 

atteint la celle du son (M = 1) 

dS/dV est positif : pour que la 


Si V» a ou M 5 1 (écoulement supersonique) =} [rte V croisse, la section doit 
augmenter 
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— On appelle tuyére supersonique, un canal qui permet d'obtenir un écoulement 
supersonique à partir d'un gaz au repos. Cette tuyere doit étre d'abord 
convergente puis divergente. Sa section minimale, ou la vitesse de l'écoulement 
> est égale à la vitesse du son, est appelée col de la tuyère. 


Onde de choc 


V<a 
M « 1 


D 4 
1 ' 
1 mmm) 
1 = Ë 


V. = do T. 
So po Pe 
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> Le gaz au repos dans les conditions pi, P T; se détend adiabatiquement dans Іа 


tuyëre. Sa vitesse atteint la vitesse du son (M = f) au col, oü les caractéristiques 
sont o, Pe, Ta V. = а,. L'écoulement continue ensuite à s'accélérer dans la partie 


. divergente de la tuyére jusqu'à atteindre le nombre de Mach M désiré ; les 


caractéristiques devenant alors p, p, T, V, a. 
.. En général, l'écoulement redevient ensuite subsonique en formant une onde de 
` choc droite qui sera étudiée dans la suite. 


Onde de choc 


A =a, Т, 
Sc ро Pc 
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3.2.2. Étude des grandeurs caractéristiques au col d'une tuyére supersonique 


En écrivant l'équation d'état du gaz parfait, la relation d'adiabaticité et l'équation 


de Saint-Venant pour les états du gaz au repos et au col de la tuyére, on obtient : 
1 T 


T, 4 2 p, (2 py" p | 2 ri 
T. а; y+Í p, y +1 D y +1 


. Dans le cas de l'air (у = 1,40), on obtient : 


1. 2.0833 ` Pe =0,634 Pe 20,528 
T; Pi Pi 
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3.2.3. Grandeurs caractéristiques dans la tuyére ou le nombre de Mach est M 


_ On écrit la relation de Saint-Venant entre le gaz au repos et l'écoulement ayant 


>» atteint Mach M, en tient compte des relations suivantes : 
V-Ma et a° -уР = y rT 
p 


`> Ainsi, on obtient : 
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> 3.3. Étude de l'onde de choc droite 


L'onde de choc est une zone de l'écoulement où l'évolution des grandeurs 
caractéristiques est très rapide : son épaisseur est très faible. Dans le cas d'une 


. onde de choc droite, on assimile le plus souvent celle-ci à une surface de 


discontinuité à travers laquelle un écoulement supersonique redevient 
subsonique. 
L'onde de choc est le siége de phénoménes irréversibles, et on peut montrer que 


— l'entropie du fluide augmente lors de la traversée de l'onde de choc. En effet, à la 


traversée de l'onde de choc, l'écoulement reste adiabatique, mais non 
isentropique (adiabatique irréversible). 
Pour l'onde de choc, on peut toujours écrire l'équation de Saint-Venant (qui traduit 


— la conservation de l'énergie dans l'écoulement), mais du fait de la non réversibilité 


du phénomène, la loi : p p” = Cste n'existe plus 


L'étude de l'onde de choc s'effectue donc en écrivant les équations de 
conservation de la masse, de la quantité de mouvement, et de conservation de 


`» l'énergie à travers le choc. 


On se limite à l'étude de l'onde de choc droite, c'est-à-dire au cas oü la surface de 
discontinuité est un plan orthogonal а la direction de l'écoulement. On герёгега 


*  parl'indice 1 les grandeurs avant le choc, et par l'indice 2 celles aprés le choc. 
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3.3.1. Conservation de la masse 
dm = p, V, dS = p, V, dS 


Volume de contróle 


On ass le théoréme d'Euler sur le v) de contróle : 


ч (p, - p») dS х = [v aas = (o, V? = p, V "Jas x шшр + p, V = p, + p, УД 


3.3.3. Poen de l'énergie 


L'équation de Saint-Venant qui traduit la conservation de l'énergie de l'unité de 
masse de fluide reste vérifiée à travers l'onde de choc : 
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3.3.4. Relation de Prandtl 
L'équation d'Euler s'écrit : 


H H 
D sS ET em WoW ш еса 
oh PV 
Or, la vitesse du son est donnée par l'expression : а? = DP 
| а а; 
. . Donc, l'équation d'Euler devient: у (v, -V, LZ EE 
oh PV V, V, 
L'équation de Saint-Venant devient egalement ` = 
2 2 
‚ A Pr re ш> = 4 a; m а; 
2 y-lp, 2 у-1р, = A еш 
ct 


kengen. s 
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~ D’où l'équation d'Euler : 


2 2 

D; p; Ур YP: a 0 

йу = Pr et ушу ы ZS S 
С ok BV, p, Vj p; V; V V 


) y (V. -V, )= 


y-1 
ere) npn dd (v, -v, Jv, +V, V, CIE =V | 


V, V, у \2]у V УУ, 2 


2 
y bah Ë d 2 my уу-у: ed J le 
V, V, 2 V, V, Y+1 2 y+1)\ 2 
kel 


1 
AË wi 


2 
Le produit des vitesses de l'écoulement en amont et en aval du choc, est égal au 
~ carré de la vitesse du son au col de la tuyère : 
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3. | 
| 
ы 22222 f enceinte dans ` 
d | E laquelle débite : 
Pi | 7 | la tuyére à la 
Pi | 
I ` | 7 


2277222222 | pression p 


Onde de recompression | 
м. | со sortie | 
Premier cas: =  )» )» ) ))))»)» DT 


V, <a, [[ж====== е Écoulement compressible analogue à celui des fluides incompressibles 
(oü a, vitesse du son au col) Ï 


La pression augmente et la vitesse diminue le long du divergent 


Remarque : La discontinuité de la pression à la sortie provoque une onde de 
recompression qui pourrait modifier les conditions de l'écoulement 
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| 
* 2222222 | E enceinte dans x 
Í | Р laquelle débite | 
A p ; p Іа tuyére à la 
^ idi | LE | pression p ` 
- | onde de recompression 1 
col sortie 


Second саз: 
V, =a, Imm» D'après le théorème de HUGONIOT, V, ne peut être supérieure à а; au 
(où a, vitesse du son au co) niveau du col de la tuyère où la section est minimale. 
Dans ce cas, la pression continue à augmenter et la vitesse diminue 
tout au long du divergent, et 3 cas sont possibles : 


(0-1) 
р2= р md | 1 peut être calculé à partir de Saint-Venant 2 (t < ) M? (a J +) 


р,>р md le jet se dilate à la sortie, et crée une onde de dépression, qui ne peut 
S remonter le courant et modifier l'écoulement 


P><p > || se produit une onde de recompression à la sotie, qui remonte le courant 
dans le divergent et provoque une élévation de la température. 
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Onde de choc stationnaire dans le divergent : 
Écoulement supersonique entre le col et 
l'onde de choc, et subsonique ensuite (c) et (d)) 
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Lx x 
V,= 0 177772 => enceinte dans - 
Т, = 3 P mede 
i ' i | 
=> V р Seege, 
Pi j _ "> la tuyère à la ` 
— T | 

Pi ” Col ЖШ — pression p 

| LII 
4 $ 


И mm) ШЕЕ 


3 
| 
> 
ta 
"du 
| (s 
M 
> |S 
EN 
И т Ny 
S |S 
— e 
SS 
M 
| 
таа 
= |s 
Wa Z fh 
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Il 
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ai 
= 
~ 
< Q` ко 
MES 
= 
— 7 
P >> TUS 
=S ^s 
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2 


EI EEN Duc Ë Jp "= Lët 


` Le débit т est maximal si X est maximal, et X est maximal si : dX/dp = 0 


WC Ren expression de la pression critique p, 
| m 
fluide incompressible 
= 2l», - p) =) - | 24 Pi 
` Donc: | | 


fluide compressible 
ES DER Gë 
T. Pi Di 


Le débit massique maximal est obtenu quand les caractéristiques du fluide à 


l'extrémité de la tuyére sont égales à celles de l'état critique (vitesse du fluide 
égale à la vitesse du son : M=1). 


Chapitre 2 


Fluide idéal isovolume : Ecoulements rotationnels 


4 z 
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FLUIDE IDÉAL ET ISOVOLUME : 
ÉCOULEMENTS ROTATIONNELS 


1. Introduction j 
En mecanique des fluides, et dans la plupart des situations, il est 


intéressant d'interpréter les événements dans un écoulement en terme de 
rotation. 
Dans ce chapitre, on se limite à l'étude de la dynamique de cet 
écoulement tout en s'intéressant aux interactions entre ce dernier et les 
| événements qui s'y associent. L'objectif, est de déterminer une certaine 
distribution de cette rotation. 
Ainsi, parler de la rotation du fluide revient à parler du rotationnel du 
. vecteur-vitesse, ou encore du rotationnel tout court. Et parler du 
rotationnel, revient à évoquer la rotation d'une petite ou grande partie de 
l'écoulement. Cette rotation se traduit par des tourbillons ou structures 
| tourbillonnaires, appelés vorticité (provenant du mot latin vortex). 
Donc, la vorticité est une formulation mathématique de la dynamique des 
fluides relié à la quantité de vitesse angulaire ou de rotation que subit un 
fluide localement. 
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Ў 2. Définition du vecteur tourbillon 


_ Soit U =(u,,u,,u; ) un vecteur vitesse. Le rotationnel de U est défini par : 


Ё и, e e 
rot (U )=VAU=| = |^ (u, ele 
Ох, Ох, Ox, 
u 
0 Í Ou, Ou, 
| Ox, Ox, OX, 
` Un écoulement est dit rotationnel, cela signifie que : 


. . Le vecteur taux de rotation d'une particule fluide, est noté o. Celui-ci est donc le 
. demi-rotationnel : 
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.. Si on désigne par [7 le vecteur vitesse d'un point M d'une région d'un fluide où 


l'écoulement se fait sans déformation (c'est à dire comme un solide), on peut 


. alors écrire : 


Ü(M)- Dot, OM 


avec: > O : Origine 
> © :taux de rotation ou vecteur instantanné de rotation 


Dans le cas d'une rotation pure, nous avons : U(0)-0 et U(M)=@A0M 


+ Si le vecteur taux de rotation œ est dirigé suivant l'axe Ox; , donc: 


0| [x - OX, EN 0 
U(M)=@AOM=| 0 |a| x, |=| ex, | et rot(U )=VAU=| 0 |=2.@ 
@) |X 0 20 


2. l1—u-.lzg > 
D'où : пе To U nr ^U  estappelé vecteur tourbillon 
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ч 3. Flux du rotationnel du vecteur-vitesse 


D'aprés le théoréme de Green, le flux du rotationnel à travers une surface fermée 


(F) peut s'écrire : 


[ror (U )л.а5 = | div Lat ао 
' F D 
Or, la divergence du rotationnel est nulle en tout point oü le vecteur-vitesse a des 
dérivées réguliëres : 
= (= =\ O (Ou ol ô (ôu, ol д (ôu, ди 
div (rot ( Ü ))= - ——— REL —-— | —-— |=0 
Ox, (Ox, Ox,] Ox, | Ox, Ox J 0х, | Ox, Ox, 
Il en résulte alors que le flux du rotationnel à travers la surface fermée (F) est donc 


nul. 
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> Remarque : 


En France, on appelle souvent œ vecteur tourbillon, mais il y a d'autres usages 


… du mot < tourbillon > ; ce qui entraine parfois des confusions. 


Dans ce qui suit, nous allons considérer un domaine (D) limité par une surface 


`> fermée (EL régulière par morceaux. Ainsi, le vecteur-vitesse est considéré dans le 


domaine comme une fonction régulière des coordonnées x; (oü i =1,2,3). 
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> 4. Circulation du vecteur-vitesse 


On appelle circulation Гаи vecteur [7 le long d'une courbe (C), dont l'élément 
d'arc est // , l'intégrale curviligne du produit scalaire [7 Ai: 


De JU di = [бат ой 7 estun vecteur unitaire tangent а la courbe (C). 
C C 


_… Le théorème de Stokes exprime que, si le vecteur J est une fonction de point qui 


est régulière dans le domaine (D), la circulation du vecteur (J le long d'une courbe 


. fermée (C), réguliére par morceaux et située dans le domaine (D), est égale au flux 


du rotationnel du vecteur (j à travers toute surface ouverte (S) située dans le 


domaine (D) et qui a pour frontière la courbe (C) : 


| Г аё@т = | | rot (Ü )i.dS = 2. | | @.ni.dS 
C S S 
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— Conséquence 1 : 
Si l'écoulement est irrotationnel en tout point d'une surface (S) ouverte et 
s'appuyant sur le contour (C) fermé, la circulation du vecteur vitesse le long 
de la courbe (C) est nulle : 


| Ü dii -0 


C 
~ Conséquence 2 : 
Si l'écoulement est rotationnel en tout point d'une surface (S) ouverte et 
s'appuyant sur le contour (C) fermé, on montre que : 


ZEN JEN 
SH | en ЗАП £dt =| SE emo TURU ue road — | |.T.dl 
"| Ot dU 2 


dt ` 
(р? (р? 
| grad | => |T.dl = | | rot | grad | — | .n.dS = 0 
C 2 S 2 


car : 


63 
FLUIDE IDÉAL ET ISOVOLUME : 


ÉCOULEMENTS ROTATIONNELS 


> l'accélération dérive d'un potentiel : 


EE 
dt 7 Ot 


=> "m PE U? _ 
-| PE Lm DAD ER — zat - 5. [Ü zat 
^| Ot 2 C 


- | 7.7 dt 
Ainsi, la dérivée matérielle de la circulation de [/ est égale à la circulation de 
. Paccélération. Si le fluide est idéal, isovolume et soumis à des forces dérivant 


d'un potentiel, on écrit alors (selon Euler) : 


— 


ply - f n grad( р ) = 0 ой f = — grad(N ) est l’accélération de la pesanteur. 


І 1 — > — = — — => 
Mais, comme : — grad( p ) = grad Z Donc: ў + grad(p/ p) grad(N)- 0 
p p 


ou encore : y —— grad( p/p+N ) Donc, l'accélération dérive d'un potentiel. 
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> 5. Potentiel des vitesses 
_ L'écoulement est dit à potentiel des vitesses dans un domaine (D) si, en tout point 
М(х,,Х,Хз) de ce domaine et à tout instant t, on peut définir une fonction scalaire 
> Q(X,,X5,X,,t) des coordonnées х;, х›, x; et du temps t telle que, sauf éventuellement 
= en certains points ou en certaines lignes ou en certaines surfaces isolées, le 
> vecteur vitesse au point M et à l'instant t est le gradient de la fonction @. 


PE grad(®) 

ob 

u, = — 

Ох, 
. . eo ob 
Ceci se traduit par : jp == ouencore: U, = —— 
Ox, Ох. 

ОФ 

u, = — 

OX, 


La fonction @ est la fonction potentiel. 
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. Le mot potentiel s'emploie aussi pour d'autres champs de vecteurs (pesanteur, par 
exemple). Dans le cas d'un écoulement avec potentiel des vitesses, à un instant 
donné, le vecteur-vitesse est, en tout point, perpendiculaire à la surface 
~ équipotentielle : 


Ф(х, ust. )= Cte 
q qui passe en ce point. 


Par suite, les lignes de courant sont orthogonales aux surfaces équipotentielles. 


— Le sens des lignes de courant est celui des @ croissants. 


Dans un écoulement avec potentiel des vitesses, en coordonnées cylindriques 
| (50,2), les composantes de la vitesse sont u, и, еї us. Leurs expresssions en 
_ fonction du potentiel s'écrivent : 


ОФ 1 ОФ ОФ 
= Ü = —. => и: = — 


' 0r r 00 Ôx, 


u 
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ч 6. Potentiel des vitesses dans le cas d'un fluide isovolume 
Pour un fluide isovolume, l'équation de la conservation de la masse s'écrit : 


div (Ü )= Zh J e 


x Të OX; 


S'il existe un potentiel des vitesses, l'équation de la conservation de la masse se 
~ traduit par : 


ОФ Oo Fo 
AD = L= = 
Or = Ôx, Or 
Dans un écoulement de fluide incompressible (isovolume), pour lequel existe un 
potentiel des vitesses Ï, le Laplacien AD est nul. 
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> 7. Ecoulement plan d’un fluide isovolume 
Puisque l'écoulement est plan (bidimensionnel), on va définir, à chaque instant t, 
dans un plan paralléle à l'ensemble des vecteurs-vitesses. Ce plan sera pris 
~ comme plan dans (x;,x;). 
Dans ce plan, la vitesse Ü =(u,,u, 0 )=(u,,u, ) est une fonction de point qui est 
réguliere. 
— Si le fluide est isovolume (qu'il soit stationnaire ou non), l'équation de la 
conservation de la masse s'écrit : 


R 


‚ (=\ Ou, Ou, Ou 
div (0]=< 4 22 Ls _ =0 
P ‚ P Ox, OX OX, Ох, Or ` ! : 
Cette équation exprime la condition pour qu'il existe une fonction w(X,,Xx»,t), dite 
fonction de courant et définie, à une constante additive près, par les expressions : 
SR 
Ох, í 
aR 
(X. 


u; 
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dy — Y .dx, + 24 dx, = —u,.dx, +и,.ах, 


X X 


suffisante pour que l'expression précédente soit une différentielle totale exacte 
~ d'une fonction de x, et de x.. 


= ò ò 
Ел би ди ү af) EE 


Ox, Ox, Or Or, ER Ox, i Ох, | Ox, i Ox, Ox, 


Ох, 


1 
De la relation de dy, puisqu'on a une différencielle totale excte (dy=0), on peut 
obtenir : 


dy =0=-u,.dx, +u,.dx, = —u, Б 2 +u, ad =й, Eu tu, dA =0 
Ox, Ox, Ох, Ох, 
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=> Ainsi, à un instant donné t, le vecteur-vitesse (u,uj, tangent à la ligne de 


courant, est normal au vecteur E , et lui-même normal à la ligne de 
Or, Ox, 


_ courant : vx, x15) Cte 


_ Ceci revient à dire que les lignes de courant ont pour équation : 


b vx, Xl; )= Cte 
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`= Calcul du débit volumique : 


Le débit volumique Q, par unité de longueur dans un tube de courant limité d'une 
= part par les surfaces cylindriques perpendiculaires au plan (x,,x;), ayant pour 
directrices les lignes de courant y, et w, et d'autre part par deux plans 
d'écoulements distants de l'unité de longueur, est donné par : 


Q, = [U 3i dt 
C 
dr 
où: =| U| et dé = bi 
dx dy 
de 


_… Donc: 8 
o = [| - lat = | Ü gd 
C 


T1 
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ÉCOULEMENTS ROTATIONNELS 


_ En coordonnées cylindriques, les expressions des composantes u, еї u, de la 
> vitesse U en fonction de v s'obtiennent facilement en remarquant que : 


21 NES 


X» X, 


u, =u,.cos 0+u,.sin 0 = 


ид РГТ Die . Sin Ө — —— 
Ох, 
` cos Ө = _ E 2 
Mais puisque : r 00 
1 Өх, Ох,» 


Donc: 
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— 8. Ecoulement 


| plan d’un fluide isovolume avec potentiel des vitesses 


Pour un écoulement plan, l'équation de la conservation de la mase s'écrit : 


Ou, ди 
ui T d ои _. 0 
OX, 0%, 
Cela implique qu'il existe une fonction de courant vx, xut) telle que : 
д 
posu 
Ox, 
а 2 
Ох, 
. Le caractère irrotationnel de l'écoulement se traduit par l'équation : 
І ôu, Ou, 
p ] Ox, ` Ôx, 
1 Ceci entraine l'existence d'une fonction potentiel dx, E telle que : 
ч ОФ 
U, = 
Ох, 
ОФ 
U, = 


i Ох, 


73 
FLUIDE IDÉAL ET ISOVOLUME : 


ÉCOULEMENTS ROTATIONNELS 


> En remplaçant u, et u, par leurs expressions respectives en fonction de Ф, on 
_ obtient : 


Co odo 
= + — = 
Or ó 


_ A partir du caractère irrotationnel de du mouvement se traduisant par ж. = 


0 


on obtient : Xr SS 


— Remarque : 


Si la condition d'irrotationnalité n'est pas satisfaite (écoulement rotationnel), il 


n'existe pas alors de potentiel des vitesses, mais il existera toujours une fonction 
~ de courant y, telle que : 
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> Compte tenu des expressions de u, еї u, en fonction de H , puis de y, déjà 
_ exposées auparavant et qui sont ré-écrites ci-dessous, on obtient : 


2 ОФ 
u; == с u, m 
Ox; et c 
n ду ns ОФ 
2 дх, Ох, 
. On peut obtenir alors les relations suivantes : 
Oy ОФ 
И! = = =Á J!" > E.O 
Ox, Or 
Oy ОФ 
i ~ u, = — = 
| Or, OX, 


> Soit, finalement les relations entre de Ф, puis de y, qui s'écrivent : 
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b 2. Ecoulement plan stationnaire d’un fluide compressible 
Dans le cas où le fluide n'est pas isovolume (compressible), l'équation de 
conservation de la masse s'écrit pour un écoulement stationnaire sous la forme : 


Ох, Ох, 
Cette relation permet de définir la fonction y = f(x,,x;), telle que : 
CL D 
Ох, 
Oy' 
= О.и 
ER PU; 


`` En fáisant le dérivées partielles des expressions précédentes respectivement par 
rapport à x, ensuite par rapport à x,, on obtient : 
д [ow] _ д?" 8 ô*(- pu, ) 
2 Ох, J- Ox, OX, i Ox, Or, 
ô (aw'\ cw (pu) 
É | 


E Ox, .Ох, x Ox, .Ох, 


Ох, 
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— Pour que, dans le cas du fluide non isovolume, la fonction de courant у soit 
mesurée avec les mémes unités que dans le fluide isovolume (et aussi avec les 
mémes unités que le potentiel des vitesses, s'il existe), on préfere définir une 

> fonction v(x;,x;, telle que : 


Cette relation permet de définir la fonction џ = f(x;,x;), telle que : 


où p est la masse volumique du fluide dans des 
conditions bien définies de température et de pression. 


Compte tenu de la relation suivante prevenant de dy = 0: 


dy =u | 22 | EA -0 
Ox, }, Ox, J. 


On peut écrire alors : 
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. Le calcul du débit volumique élémentaire par unité de hauteur dQ entre deux 
ignes de courant de côtes w et v,*dy donne : 


dQ, =u,.dx, —u, aray --P [88 Ea- Po. 
p V Ox, Ох, 


. est finie, la valeur de p varie le long de la courbe (C) Ag 
qui joint les 2 points M et N (voir figure ci-contre). 
Dans ce cas, le débit massique par unité de longueur 
_ (ou de hauteur) Q à travers la courbe (C) est : 


| Q = |40 = | p.40, - | oU ñ.dt - -| po dy 


= od; -v,)- p.v, р), -v; o 
Ainsi, pour un fluide à masse volumique variable, la cotation des lignes de courant 
{ y ne permet que le calcul du débit massique, et non pas du débit volumique. 
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— Dans le cas d'un écoulement avec potentiel des vitesses, les lignes de courant et 


les lignes équipotentielles forment toujours deux familles de courbes 
orthogonales. 


— Considérons une maille élémentaire constituée par deux lignes de courant 


voisines y, еї y;+dy et deux équipontentielles voisines Ф, et Ф,+аФ. 


Ob p, Ov 


> La relation u, =— =>. =Á D, +аф w, һау 


| Ox, p Ox, 
donne, en remplacant les axes Ox et Oy 
. respectivement par Pret Pn , déjà définis ci-contre : 


OD _ дү 

ou bien Or р дп 
LÀ 
PU; АУ ЖЕ» 
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> 10. Ecoulement de révolution 

Par analogie avec ce qui a été fait pour l'écoulement plan, nous définissons la 
fonction de courant dans un écoulement de révolution à partir de la loi de 
~ conservation de la masse dans un tel écoulement, exprimée en coordonnées 


cylindriques. 
La divergence du vecteur p.U s'écrit en coordonnées cylindriques : 


(00): ара ) Aen, | ‚ дри) 


Ох, 


Or 00 


> Par suite, l'équation de conservation de la masse 


n + div(p.U )- 0 


> s'écrit en coordonnées cylindriques : 


2 p1 k du), Aen) 


ôt rV дг 00 ду, 
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— En écoulement stationnaire de révolution oü, par définition, les dérivées par 


rapport à Ө sont nulles, l'équation de de conservation de la masse en 
coordonnées cylindriques, se réduit à : 


. Dans cette équation, la vitesse circonférentielle и, n'y figure pas, et le mouvement 
dans un plan méridien est indépendant du mouvement circonférntiel. Ainsi, 
l'équation de de conservation de la masse dans un écoulement stationnaire 
. méridien (u, = 0) est donc celle qui est ci-dessus. 


Cette équation de conservation permet de définir une fonction de courant y(r,X3), 
telle que : 
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… Les constantes ou еї r, sont introduites uniquement pour des raisons 
d'homogénéité. Dans le cas d'un fluide isovolume, les équations précédentes se 
réduisent à : 


. Ces équations sont aussi valables pour un écoulement instationnaire, et la 
fonction de courant v(r,x,,t) est alors dépendante du temps t. 


` Ainsi, l'équation v(r,x,,t) définit la méridienne d'une surface de révolution qui est 
une surface de courant. Si l'écoulement est méridien, cette méridienne est alors 
une ligne de courant. 
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`> Les équations suivantes : 


M AV x V P 
= UR = 
> donnent : Ur Poto Her Jom 
nn = 0 
r OX; 


`> Cette relation exprime que, dans un plan méridien, la vitesse U = (u,,x;), tangente à 


la ligne de courant, est normale au vecteur | ôv ду | lui-même normal à la courbe 
Cte òr Gr, 
: y = А 


Si l'écoulement méridien admet un potentiel des vitesses @(г,2,1), les composantes 
_ de vitesse u, et us sont données par : 
E oo 
` Or Ох, 


u 
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~ Compte tenu des expression précédentes de u, еї de u, dans le cas d'écoulement 
= stationnaire et d'un fluide non isovolume, on а: 


EP 2 = Ж. 


. 3 . 
Or Polo py. OX; 
OW рт „ 2 PT ob 


Ox, On рњ ôr 


> Dans le cas d'écoulement stationnaire ou non et d'un fluide isovolume, les 
équations précédentes se réduisent à : 


he CL 
i Or r, ôx, 
Oy r ОФ 
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. Compte tenu des expression précédentes de u, et de u, dans le cas d'écoulement 


> stationnaire et d'un fluide non isovolume, on obtient : 


pu A 2 
Or r, Ox, 
Oy т ОФ 


Ox, 1, Or 
En particulier, si en un point, l'axe des r est tangent 
à l'écoulement (axe des т) , et si l'axe des x; 
— coincide avec la normale de cette ligne (axe des n), 
on a dans le cas précédent : 
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Dans le cas particulier d'un écoulement avec potentiel des vitesses d'un fluide 


> isovolume, l'équation de la conservation de la masse s'écrit : 


EE 


r 


ôr r 00 Ox 


ОФ OD ,, , - i 
> Compte tenu deu = — et de u; = SCH l'équation précédente devient : 
r X3 


ob 1 o 0Ф a 


д =r ôr ër 


. C’est le Laplacien de la fonction Ф en coordonnées cylindriques. 
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— Dans le méme cas d'écoulement, d'apres : 


p < oo 
< и, = —— 
Or 
ob 
Из = —— 
Ох, 
… la condition d'irrotationnalité ôu, = Ou; —() permet d'écrire : 
| Or 
| 


I> Le premier membre n'est pas le 
Laplacien de la fonction w en 
coordonnées cylindriques. Le 


= | ep résultat est donc tout à fait différent 
` compte tenu de : Ог г Ox, de celui auquel conduit l'emploi des 
Я Oy т Ф coordonnées  cartésiennes “dans 


Ôx, h 6r l'étude de l'écoulement plan. 


Chapitre 3 


Représentation des écoulements plans, isovolumes et 
irrotationnels par des fonctions analytiques 


4 Ў 


| Chapitre III 87 
REPRESENTATION DES ECOULEMENTS 


| PLANS, ISOVOMES ET IRROTATIONNELS PAR 


DES FONCTIONS ANALYTIQUES 


> 1. Introduction 


On dit que l'écoulement est irrotationnel lorsque les particules fluides ne 
subissent pas de rotations pures. Autrement dit, le vecteur tourbillon О 


1 est nul dans un écoulement irrotationel. 
z Ja = Á 
О = 2 VAU = 0 
D'un point de vue mathématique, la relation : 
| VAU = У A grad(D)=0 v Ф 
est toujours vraie. 


. Onestalors en droit de définir un scalaire Ф tel que : Ü = grad(®) 
où @ correspond donc à un potentiel scalaire : c'est le potentiel des 


vitesses. 


REPRÉSENTATION DES ÉCOULEMENTS - 


. PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


- Il est donc possible d'exprimer les composantes du vecteur vitesse à 


> partir du potentiel des viteses : 


ОФ 

U, = —— 

OX, 

__—_ ОФ 
U = grad(®) => 4и, = — 
OX, 

ОФ 

Из = — 

OX, 

~ L'écoulement étant incompressible, on doit alors vérifier : 


div (U)- 0 
. ce qui conduit à la relation de Laplace : 
ôP Go 0° 


EE ou A(D)=0 
Ôx; Ox; дк @) 
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^ PLANS, ISOVOMES ET IRROTATIONNELS PAR 
| DES FONCTIONS ANALYTIQUES 


ь 2. Définition de quelques propriétés des fonctions analytiques 


Une variable complexe z = x; + ix, peut être représentée dans le plan 
(XX), par un point Mz(x,,x;). z est appelé l'affixe du point M, ou tout 


simplement le point z. 


. La distance OM =r=.Jx; +x est le module de z. 
. L'angle @que fait OM avec Ox, est l'argument de z. 


. Donc : 
| š ES d 
esa sp. = r(cos O+isin 0)e 


| avec: 


REPRÉSENTATION DES ÉCOULEMENTS " 


PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


— Dans un domaine (D) du plan des х;х,, toute fonction f(z) de la variable 


complexe z peut étre écrite sous la forme : 


f( Le 9, ox, )+iwlx,, x, ) 


. qui met en évidence la partie réelle Ф et la partie imaginaire w. 


Ainsi, une fonction f(z) est dite analytique au un point z, si l'on peut 


trouver un cercle de centre z et de rayon non nul tel que, à l'intérieur de 


~ се cercle, les fonctions @(x;,x;) et y{x;,x.) soient univoques et aient des 


dérivées partielles continues vérifiant indentiquement les relations, dites 


. de Cauchy-Rienmann : 


REPRÉSENTATION DES ÉCOULEMENTS е 


à PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


— La différentielle de f(z) s'écrit : 


df = E Jas, + (2) as, (2и Jas, (EM 
x, Ce x; x 


Ob OV et? 


> Les relations `". ОСИ permettent d'écrire : 


Ox, Ox, Ox, x, 

ob 0 
df =| E «i-o. ld 

Ox, Ох, 

Ainsi, quelque soit la variable z, on obtient la méme limite : 

d 
df € oo + l дү 
dz Ox, Or 


L'existence et la continuité des dérivées d'un ordre quelconque de f(z) 


` permettent de déduire, à partir des relations ci-dessus, les équations : 
| dO o ду © 

-4 =0 et Ee 

Od ` Ôx? Óx; WS 


REPRÉSENTATION DES ÉCOULEMENTS - 


| PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


> 3. Potentiel complexe et vitesse complexe 
1 Soit, au point /М-(х,,х,), un écoulement plan à la fois isovolume et 


irrotationnel. A chaque instant, les composantes u, et u, de la vitesse 


- О(х;,х,) en un point, satisfont les relations : 
ОФ _ дү ОФ ` _ дү 


U, = et и, = 


ex, m Әх, Ох, 
La fonction complexe f(z), définie à chaque instant par l'expresion : 


f( z) = dX x,, x, )t iv( x,, x, ) 
est au point z, une fonction analytique. Nous l'appellerons potentiel 
_ complexe au point z. 
`La dérivée de df/dz, qui est aussi une fonction analytique au point z, est 
donnée par la relation : ob ð | 
нба) а tia u, — iu, 

< X 

Cette dérivée, appelée vitesse complexe w au point z, est l'imaginaire 
conjugué de la fonction complexe U = u,+iu, représentant le vecteur 
vitesse au point z. 


REPRÉSENTATION DES ÉCOULEMENTS °° 


é PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


a) Ecoulement uniforme 


Considérons l'écoulement plan modélisé par le potentiel complexe 


des vitesses : f(z)-uz 


. 
|] | Onaalors: Ф(х,у) +/ W(x,y)= Ш(х +iy) = Ох + iUy 
| X = Ux 
Par identification, on obtient : 5 У) 
: 4 Les lignes de courant sont telles que : Y{(x, y) = Uy = C* 
> y =C" Vx ce sont des droites horizontales. 
l Les équipotentielles sont telles que : ф(х, у) = Ux CTT 


> Xx —C' уу ce sont des droites verticales. 


Détermination du champ de vitesses : V = 


: : - = У 
La vitesse est uniforme : v = Ue, oy ox 


Ay HY A A AM | 


REPRÉSENTATION DES ÉCOULEMENTS o 


Í PLANS, ISOVOMES ET IRROTATIONNELS PAR 


DES FONCTIONS ANALYTIQUES 


lignes de courant : V(x,y) = Uy =C" = у =C" vx (droites horizontales) 


H équipotentielles : g(x, y) = Ux = C^ > x = C* vy (droites verticales) 


champ de vitesses : V = Ue, 


g = С° 
w = Cte 


écoulement uniforme 
f(z) = UZ 


Ay A HY AN AM С 


REPRÉSENTATION DES ÉCOULEMENTS е 


é PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


b) Ecoulement plan autour d'une source ou d'un puits 


Considérons l'écoulement plan modélisé par le potentiel complexe 


des vitesses : | . D 


et Cune constante réelle. 
=> f(z) =Cln(r e'*) = Cnr +i6) 
On peut alors en déduire la fonction de courant et le potentiel des 
vitesses : He PESE 
Kafe AER e, 
Les lignes de courant sont telles que ` V(r,80) = СӨ = C": 
=>  0-—C' vr ce sont des droites passant par l'origine. 
Les équipotentielles sont telles que : @(r,0) = CInr = CT 


— 


=> г=С©® Y8 ce sont des cercles concentriques 
centrés sur l'origine. 


REPRÉSENTATION DES ÉCOULEMENTS °° 


PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


Détermination du champ de vitesses : 


E vica 9" n F: 
_10p _ OT 
^ roO ar 

. (v = С/г 2 8 

Soit : V = |» x sy <e 


La vitesse est donc radiale et 
T inversement proportionnelle 
ф =C à la distance à l'origine. 


Si C>0, alors l'écoulement est dirigé vers l'extérieur 
=> écoulement divergent => source à l'origine. 


a Si C«O, alors l'écoulement est dirigé vers l'origine 
=> écoulement convergent = puits à l'origine. 


AE HY A HY HY | 
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PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


Signification physique de la constante C : 


Calculons le débit volumique de cet écoulement radial (source 
ou puits) : 
MET H V.ndS оо S est une surface fermée 
= entourant l'origine. 


L'écoulement ayant lieu dans un plan 1 à l'axe z, on peut considérer 
comme surface d'intégration un cylindre de hauteur Az=1, et donc : 


H Ae = Asa Il reste alors à intégrer sur un cercle de 

: rayon r quelconque, centré sur l'origine. 

" r. — V = С/ге 

q, = Az фӯлгаө = azr | V.ndó ой É Si ў 
£ 0 n-e, 


2; 2; Se . 
> q, -azr| "Сав лг | "de =2zCàz débit volumique par 
o e r Jo unité de hauteur 


Inz 9,70: débit de la source 


q А q 
>L- wv = = v 
7 et par conséquent ` f(Z) = cn ER s 


Ay HY HY HY AM С 
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é PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


c) Vortex ou tourbillon libre 


Considérons l'écoulement plan modélisé par le potentiel complexe 


des vitesses : | E "m 
ARE où z-x+iy-re 
et C une constante réelle. 
=> f(z) --iCin(re?^)- -iCc(nr +i8)=C6-iCinr 
On peut alors en déduire la fonction de courant et le potentiel des 
vitesses : e(r,80) = Ce 
V(r,0) = -Clnr 


Les lignes de courant sont telles que : Y(r,0) = —Clnr = C* 


E te š 
=>  r-—C" VÓ ce sont des cercles concentriques 
centrés sur l'origine. 


Les équipotentielles sont telles que ` Ф(г, 0) = СӨ - C* 


= Ө=С©® Vr ce sont des droites passant par l'origine. 


REPRÉSENTATION DES ÉCOULEMENTS `` 


PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


Détermination du champ de vitesses : 


| : | et inversement proportionnelle 
ECS à la distance à l'origine. 


Si C>0, alors l'écoulement s'effectue autour de l'origine dans 
le sens trigonométrique. 


Si C<0, alors l'écoulement s'effectue autour de l'origine dans 
le sens horaire. 


LJ 
|] N | Я La vitesse est donc orthoradiale 
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PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


Signification physique de la constante C : 


Calculons la < circulation > de la vitesse autour de l'origine : 


T = bv. où dí parcourt une ligne de courant 
{ quelconque, Le une cercle de rayon r. 


_ C. 0 _ 2x 
Avec : V= — €, et d =rd0e, > zi “габ =2лС 
0 
Donc : E EA 
27 


Р E à oü Г est la circulation du 
et par conséquent : ©) = ~ Inz | | 
P q (z) 2л vortex (tourbillon libre) 
Si Г>0, le vortex tourne dans le sens trigonométrique. 
Si «0, le vortex tourne dans le sens horaire. 


Ay HY HY HY AM С 


( On a alors : | 


REPRÉSENTATION DES ÉCOULEMENTS `" 


Í PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


d) Coins et points d'arrét 


Р ` On appelle < point d'arrét > un point oü la vitesse est nulle. 


Considérons l'écoulement plan modélisé par le potentiel complexe 


ë des vitesses : : 
EE EE 


En coordonnées cylindriques : z - re^ et donc : f(z) = C r^? 1e"m-De 


e(r,8) = Cr"? cos (m + 1)6| 
Vir. 0) = Сг" sin [Gr 4 1)6] 


Г op Lon 

š , = FO d er I r o8 
M Le champ de vitesses s'obtenant par : V E 12ф _ _ dP 
^ гәө or 


v, = C (m + 1)r” cos [m + 1)8] 
v, = —C(m +1)r” sin [m + 1)6| 


— 


On trouve : v = | 


On remarque que v, = v; = 0 роигг= 0 = l'origine est point d'arrêt. 


REPRÉSENTATION DES ÉCOULEMENTS `” 


é PLANS, ISOVOMES ET IRROTATIONNELS PAR 
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La ligne de courant passant par le point d'arrét doit donc vérifier : 
| Ҹ(г,0) = Сё = Ф, où W, = W(r,,@A) = Cr sin[(m + 1)0, |= 0 
L'équation de cette ligne de courant s'écrit alors : 


B - point d'arrét 
! = 


nhi si = 
cr віп[(т + 1)0|= 0 = fe + 1)0]= 0 vr 


avec Пе 


sin=0 : 0-0Ovr => demi-droite Ax 
sin=1: Өө = == = а Үг => demi-droite d'angle с avec Ax 


sin-z2: 0- E 2o vr => demi-droite d'angle 2« avec Ax 


m + 1. 


. 
п | 
| » g= л Vr = (m+1)9 = nz vr 
i (m +1) 


...etc. 
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é PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


| Les lignes de courant pouvant étre 

assimilées à des barrières infranchissables, 
» celles passant par le point d'arrét forment 
des « coins » : ce sont les coins d'arrét. 


Analysons maintenant l'écoulement 
| du fluide entre ces coins d'arrét pour 
» quelques valeurs particuliéres de m. 


L] 
| |] ` e p Z EN 
` 7 TS ` 
EN A s x À . Ei 
py M \ 
` A1 
d B 
4 BU Pul Pg LE Í Í PL < 


4 
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PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


Hz) = Cz" où me 


cas oü m=1 
л = 


EIN coin à angle droit 


Ҹ(г,0) = Сг? ѕіп[20|-с'° et EE 


= Ҹ(г,0) = 2Сг? ѕіпӨсоѕ0 = 2Cr sin гсоѕ0 = C* 
Е N = 
y x 


x 
— 


tr, gi EE en 2Cxy-c* 


c'e à l'intérieur de ce coin 
У = — les lignes de courant 
sont des hyperboles 


Les équipotentielles étant en 
tout point 1 aux lignes de 


к ОООО x pi p узд également 


tt 
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PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


cas oü m» 1 = тыс cm 5 <a=_ ат 
Z 
V Z 
Z 
С 
2 
0222 х 2222222 
cas oü 1<m<0 EE cas ойт=  «-2z 


| У^ 


Z > Zx уу H ve —— sa ^x 
| 
j | 
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é PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


e) Doublet et dipóle 


On a vu que pour qu'un écoulement puisse étre décrit au moyen 
d'une fonction de courant v et d'un potentiel des vitesses o, il faut 
que ces deux fonctions vérifient l'équation de Laplace : 


AP —0 et Ap=0 > f(z)=o+iW 
Considérons alors 2 écoulements tels que : 
(1) AV, =0 et Am =0 > А(2) = ф +i, 
Comme l'équation de Laplace est linéaire : 
A(194 + 7292) = MAP + Аф) = О 
NES A + 2%) = A AU. F AAP = О 
Donc, si on pose: @ = Am + оф, et W = AV + ZY,, alors 
AY =0 et Ag =0 = f(Z)-9« iP = A,R(Z) + 2 (2) 


Par conséquent, f(z) décrit l'écoulement résultant de la superposition 
des deux écoulements f, et f. 


Ay AN HY AN AM С 
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PLANS, ISOVOMES ET IRROTATIONNELS PAR 


DES FONCTIONS ANALYTIQUES 


On peut donc superposer plusieurs écoulements élémentaires pour 
créer des écoulements plus évolués, et ceci par simple addition des 
|] { potentiels complexes correspondants. 


Association d'une source et d'un puits : 


| |] Р Considérons une source de débit +q, située еп x-a, à laquelle on 
superpose un puits de débit -д, situé en x=-a. 
Le potentiel complexe résultant s'écrit : 
| К ыа 
f(z) = + In(z - a) _ 9 [n (z + a) Posons : xL Asi Gef 
27 Z, = 2 на = пе! 
source puits 
d | рой: f(z)- Z-(Inz, -Inz,) - 9 (Ina «ie, —Inr, — i0,) 
2л 2л 


r 
SE 


et 
Ë => f(2)= Linh + i(0, Z = 2л T, 
| E Y = 9 (0, — 6) 
| 2л 
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PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


p= Tin Et donc, les lignes de courant sont telles que : 
25 d ei ` 


| Lë Ecoulement généré 
par un doublet 


Les lignes de courant sont des cercles passant tous par P et S. 
. Les équipotentielles sont également des cercles. 


REPRÉSENTATION DES ÉCOULEMENTS `” 


é PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


Faisons tendre la distance entre le puits et la source vers 0. 


ы ИШИНЕ! 
x 
| fz) = +9102 a) in + a) = Lin = Z MERA 
f(z) = 5; 9 1022/2) oü ar a ө” 
Donc f(z) = Zn - a/zY] = Z 2ln(1 - a/z) = Z x 
л 


Ay AN ÜA HY LAM С 


À Posons 2aq = p le moment dipolaire : "m 


REPRÉSENTATION DES ÉCOULEMENTS ' ^ 


Í PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


ee = _ JL P (cose-isine)- + iw 
| 2л r 


1 p 1 
p=—- cos 0 vV-C* о — P sing = c' 
D'où 27 r S 2л r 
E Y = + Bang équation d'une 
2л г 


ligne de courant 
| — Zsin@ =c" = гіпо = С‘ r? — y = C* (x? + y?) 


=> x*'4y*-C*y-0 —x^-«y*-Ky-0 > x2+(y-K/2) = (K/2Y 
been — 
Е équation d'un cercle de 
Y-C" o pesas (0,K/2) et de 
rayon K/2 
Les lignes de courant sont des cercles 
tous centrés sur l'axe y, et passant 
tous par l'origine. 


Ay AN Ü A AM ` 
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PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


Ecoulement généré 
par un dipóle 


| 
APENS x 


A Remarque : 
| Le vecteur moment dipolaire p — PS donne l'orientation globale de 
l'écoulement. " 


| REPRÉSENTATION DES ÉCOULEMENTS Wës 
PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


>» 4. Ecoulement autour d’un cylindre circulaire 
Considérons l'écoulement résultant de la superposition d'un écoulement 


uniforme, parallele à l'axe des х,, dont le potentiel complexe peut s'écrire 
~ fi(z)=U .z, et d'un doublet plan à l'origine, parallèle à l'axe des x, et dont 
le potentiel complexe f(z) est mis sous la forme : 


dé 2 
J. (z )= 


> Le potentiel complexe résultant a pour expression : 


паела. 


II peut s'écrire : 


2 2 
a 
панни и) ш 6 


Z 


= REPRÉSENTATION DES ÉCOULEMENTS We 
PLANS, ISOVOMES ET IRROTATIONNELS PAR 
0 DES FONCTIONS ANALYTIQUES 


~ ou encore sous la forme : 


ж ia .x 
JG)» iy ets +з = d 


2 2 2 
X; + X> X; + X> 


a° | a° 
zu JPE АО 
X E X T4 


>Le potentiel des vitesses et la fonction de courant sont respectivement : 
~ 2 2 2 

` a a a 
| b=U x, À I+—— |= U..r.cos Ө. =; =U -cos 0. r += 
x° + x r 

bk 2 2 

| а . a° . a 

y =U Xo) 1-——|=U,r.sin0] 1-— |ZU,.sinQ. r — 

X; tx, r r 
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PLANS, ISOVOMES ET IRROTATIONNELS PAR 
i DES FONCTIONS ANALYTIQUES 


Ë un point z, la vitesse complexe est donnée par : 

L a° a° Ф 
wz)- Ï = L Ua ener NEL E 
8 dz dz S Z Ox, Ф, 


donné par le potentiel complexe U, z. Il est donc uniforme, parallèle à 


= l'axe des х;. 
= — 
F О. 
— = 
p —+ 
—s — 


1 REPRÉSENTATION DES ÉCOULEMENTS We 
PLANS, ISOVOMES ET IRROTATIONNELS PAR 
DES FONCTIONS ANALYTIQUES 


— La ligne de courant particulière qui correspond à w=0 est constituée par 
l'axe des х,(0=0) et par le cercle de rayon a ayant son centre à l'origine. 
Pour 2=+а, c'est-à-dire aux points d'intersection A et A’ du cercle de 


. rayon a avec l'axe des x, la vitesse est nulle ; ce sont deux points d'arrêt. 


» l'infini amont au point A’ et du point `. 


© La relation w(z)-df(z/dz-U,(1-a?z?) ` _ 
montre que, sur l'axe des x, de . 


_ A à l'infini aval (/z/>a), la vitesse est 
positive, tandis qu'elle est négative 
» sur le segment А'А (/z/<a). 


. On peut donc considérer que dans le sens positif, la ligne de courant y=0 
> est constituée par l'axe des x, sauf la portion AA’ où cette ligne de 
courant se dédouble en deux-demi cercles A'BA et A'B'A. 


` A REPRÉSENTATION DES ÉCOULEMENTS "P 
" PLANS, ISOVOMES ET IRROTATIONNELS PAR 
| DES FONCTIONS ANALYTIQUES 


— Les lignes de courant situées à l'intérieur du cercle de rayon a se ferment 

_ sur elles mémes, de sorte qu'une particule fluide intérieure reste à 
l'intérieur. De méme, une particule extérieure au cerle reste toujours à 
. l'extérieur. On peut alors remplacer le cercle de rayon a, qui est ligne de 


I courant, par un obstacle solide fixe de méme contour. 


. La vitesse en un point de coordonnées cylindriques r et 0, s'obtient à 
| partir de l” expression du potentiel scalaire Ф, déterminé précédemment : 


2 
a 
P =U co det 
r 
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Á К 
4 77 
г 


- En particulier, pour un point de la circonférence (r=a), on а: 


s 


2 
a 

u -Pu col À = 0) 
a 


A 


I ^. др 
4 100 а 
4 Ио == UI. em O. 1+ = —2U .sin0 
г 00 
Ainsi, la vitesse U au point d'azimut da pour moque. 
U =2U {sin | 


7 | 
B et B’ où le cercle coupe l'axe des — 


u А H i 
X, Elle est alors égale au double — 
T Be. í 


' de la vitesse à l'infini. 


` REPRÉSENTATION DES ÉCOULEMENTS "° 
" PLANS, ISOVOMES ET IRROTATIONNELS PAR 
j DES FONCTIONS ANALYTIQUES 


> Si nous superposons à l'écoulement précédent, de potentiel complexe : 


пои) 
Z 


. un vortex placé à l'origine, et pour lequel la circulation est Г, nous 


obtenons un mouvement dont le potentiel complexe est : 
2 


патен s Tue) 
Z 2л 


. Dece potentiel, on peut tirer : 


2 2 2 
а.х Гө a | T a | T 
b=U x; += E Оса Et => EU m 0 r+ = ni 
x FX; | 27 r 2л rj 2m 


2 
a. o E r 


a x ; a 
y =U It = RN КҮТ. = NU. pz MU 
X 2л г) 2m 
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— Concernant les lignes de courant, le potentiel des vitesse augmente de Г 
pour un tour effectué dans le sens positif autour de l'obstacle. 


. Les lignes de courant sont alors symétriques autour de l'axe des х„ mais 
non symétriques par rapport à l'axe des х,. 


L'expression de la vitesse complexe : 
m k 


d a Ш 
ieu z — =u, — iu, 
| dz Z ZG 
. Montre que cette vitesse s'annule pour : 
z- 2-а? =0 
q 2zU_ ES 
Le discriminant de cette équation est: — _ +44? 


sik 


REPRÉSENTATION DES ÉCOULEMENTS Wes 
PLANS, ISOVOMES ET IRROTATIONNELS PAR 
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3⁄4 A 
Ss EP 


~ Trois cas peuvent se présenter : 


D a) Premier cas : < 4,7 


>a 
П existe deux points A et А’ de la vitesse nulle, définis par : 
il d 

+ 


2 


Z, =X, +X, = —— + |а ————— 
1 7172 дщ], 167U; 
2 r 2 r 
dup RU == > Ap = [Œ = 2 
- 16.2U 5 16." U 
r x Г 
| 2 4л] " Asti 
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>- Ces deux points ont une ordonnée commune 7/7(4zU,) et sont sur le 
_ cercle, (x,? + x,2 = ai. х, 

La ligne de courant qui passe par 
— les deux points A et A est : 


y =v, = La) 
Cette ligne de courant joue le 
. même rôle que celui de y=0 dans 
le cas circulaire traité 
… précédemment. 
> Ce n'est pas une ligne simple : elle comprend à la fois deux branches 
infinies et le cercle de rayon a. Les points A et А’ de vitesses nulles sont 
` des points doubles dont les arguments 0, et 7-0, sont données par : 


Cette expression montre la correspondance entre la 
valeur de la circulation autour du cylindre et la 
position des points d'arrêt A et A’. 


X4 
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b) Deuxième cas : = 4,7 

A ya 

Les deux points A et A de vitesse nulle, viennent se confondre sur l'axe 
1 des x.. | X; 


2 sS | | 
" 
> 
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Во) Troisième cas : 

. Les deux points de vitesse nulle, sont 
“sur l'axe des x, L'un des points est à 
"l'intérieur du cercle, hors de 
… l'écoulement que nous considérons. 
> L'autre, E, extérieur au cercle, 
. appartient à une ligne de courant qui 
| présente une boucle à l’intérieur de 
"laquelle les lignes de courant sont 
- fermées et entourent l'obstacle. La 
. boucle limite alors une < atmosphère > 


fluide autour de l'obstacle. 


Chapitre 4 


Dynamique des fluides visqueux : Pertes de charge 


4 T 
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DYNAMIQUE DES FLUIDES 
VISQUEUX : PERTES DE CHARGE 


1. Notions générales 


Étudier l'écoulement d'un fluide réel revient à résoudre l'équation de Navier- 
Stokes. 

En pratique, cette équation aux dérivées partielles, ne peut se résoudre 
analytiquement que si on poses des hypothéses simplificatrices, en distinguant 
deux grands types d'écoulement : en régime laminaire ou en régime turbulent. 


poeni [кик [еюн 
régime laminaire „ ^ régime de transition régime turbulent 


BE [Farc ans, «usan, 08] M prm esa 
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poen jon 


régime laminaire régime turbulent 


En régime laminaire, on pourra généraliser | En régime turbulent, on utilisera 
l'équation de Bernoulli en introduisant la | des relations empiriques 
notion de pertes de charge dues à lal généralement obtenues à partir 
viscosité du fluide en écoulement. de l'expérience. 


La transition d'un régime d'écoulement à un autre est caractérisée par le seul 


paramètre adimensionnel représentant le nombre de Reynolds et noté Re. Celui-ci 
est défini par : 


: masse volumique du fluide en écoulement (en kg/m?) 

: viscosité dynamique du fluide en écoulement (en Pa.s) 

: viscosité cinématique du fluide en écoulement (en m?/s) 

: vitesse débitante de l'écoulement (en m/s) 

: diamètre équivalent de la conduite (en m): diamètre hydraulique 


D < < = D 
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La transition d'un régime laminaire à un régime turbulent s'observe pour : 
Re = Re, = 2000 (nombre de Reynolds critique) 


transition Re 


On peut étendre le domaine oü le nombre de 
Reynolds Re d'un régime est laminaire au-delà de 
Re, à condition de prendre certaines précautions 
expérimentales (sans perturbations) 


v Quelles que soient les perturbations subies par le système (dispositif 
expérimental), le régime reste laminaire si Re < 2000. 


` v^ Pour 2000 < Re < 4000-5000, le régime est dit de transition. 


v Le régime turbulent est atteint à partir de Re > 4000-5000. 
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> 2. Ecoulement laminaire — Pertes de charge régulières 


` Pour un fluide newtonien incompressible, l'équation de Navier-Stokes s'écrit : 
dV E => " 
Po VISUM ande 
Si l'écoulement est permanent (stationnaire), on a : 


fz) 
Е AX + (v.v)V = (MERDE V A(V AV) 
0 


it : 

| Soit 1 
ou bien 
ou bien : 


| oubien: 
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DYNAMIQUE DES FLUIDES 
. VISQUEUX : PERTES DE CHARGE 


Zo => 1 ips => = 
1 (ptor ov =u AV +2pV nL 


Projection sur une ligne de courant : 


" 1 _ 2 
donner un - EE 


A BE 
= 


(v ^ à) e, =0 
| - 1 2 _ = _ 1V UV 
>> V p+pzg ЕДА ds е, =u AV dee, 
. Donc, la projection sur chacun des trois axes donne : 
д 1 2 Е 
ds ë, = dx, ë, + dx, ë, + dx, ё, = [psp ase lov ШЕ 


3 


AV = Au, ë, + Au, ё, + Au, ё o 1 2 Ss 
1 61 me 3 Єз Zens TL ES и Au, Ù, 


o А = Ô 3 Ô = Q _ Q | ] d 
pere EE EE E ME EE = и Au 
1 EN К j E i Ох, 2 м, Kä 
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D'où l'expression finale : 


Ô 1 
< [p+pgz+ZpV*]=p am 
Ox, 2 
д 1 
2 [р+рв e+ Aer" jean, 
Ox, 2 
Ô 1 
2 (p+os z+Z pV? ) p Au, 
| Ox; 2 
. Pour un écoulement laminaire unidirectionnel s'effectuant le long de l'axe x;, on a: 


u, =V 
_ | 1 
|. V4u,-0 => V=Vë, etsion pose: p+pgz+>pV =E 


CP. и Au 


ôx, | Pr [x,X2, Xa )= Pr (x) 
RK Op, ` 2 2 2y 
| = Ox, ° Е “ла u AV =u aR 
e » Op, = dx; ôx? ] Or x 3 
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CR D'aprés l'équation de continuité, on a: 


0 0 
= OV 
VV=0 e nee. ES en -0 < MEY e =0 


| parcourue par le fluide. 


7 = 
ü Ox, (Ox, Ф, Ox, Ox, OX; 
D'où : 
Fonction de x; Fonction de 
X; et de x; P: 


On en déduit alors que la charge 
varie linéairement avec la distance 


X 
d . = т L4 = 
Pt Cet eg n La charge diminue linéairement avec la 


dx , progression de l'écoulement de fluide. 
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dp H." 4 Ap 
E = = Í 
P а dy "x x. x "x, 


i 
ps u E 


v v M 
Pression totale en (1) Pression totale en (2) perte de charge 
régulière : Ap, > 0 


II ne reste qu'à caractériser le gradient de pression totale (ар, / ах). 
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3. Ecoulement de Poiseuille 
. Soit une conduite cylindrique, de rayon R, disposée horizontalement, et dans 
laquelle s'écoule un fluide isovolume (incompressible) en régime laminaire : 


y 


V=Vé, > У, =V,=0 


L'équation de Navier - Stokes se réduit à : 
dp 
_ = H AV 
dx, 


. D'autre part, l'équation de continuité s'écrit : 


К 0 0 
VV =0 => 106), 19€ ди, y — Op OV _, 


r Or г 20 Ox, Ox, Or 
| La symétrie de révolution de la conduite impose que : 
д OV 


—=0 = —= 
00 00 
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`< 
`< 
E I 


», D'où le Laplacien s’exprime sous la forme : 


| 

1 -12( +) V M -12( +) 
il x il x) 

aw 


fonction de x, fonction der 


À Par une simple intégration, on obtient le profil de vitesse V(r) = u,(r) : 
» BE A dE A dv Ar dV Ar B 
Е = => B 

Or dr 


H dr H dr u2 dr u2 r 


r Or 


4 
4 ^ 


(Sea 
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. Détermination des constantes B et C grâce aux conditions aux limites : 
2 
AR ‚ 
r=R = V(r=R)=0 = ШР A mmm». adhérence 


> |\r=0 > V(rz0)-Cte (finie) = В In(0)- C-Cste = B=0 

D D'où: 

| C- A К Е т 5 profil de vitesse 
" = д 4 parabolique 


EM = Pour avoir V(r) > 0 quel 
= que soit r, il faut que A« 0 
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- — Débit volumique Q, : Е "NT 
dQ =Ү(г) 45 -2nV(r)rdr > Q, = [27 V(r) as [a | (R Ch dr 
^ r=0 r=0 


"ү [ает А Ed 
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=> Vitesse moyenne de l'écoulement V, : 


_ T 2 
7 _0,_ 128 и \ dx, 1 dp, | p 
" s$ D' 32 u SE 
"4 


х, 2x; 


1 Pour une conduite de longueur L, la perte de charge totale Ap, est donnée par : 


dp dp 2) 
Ap =p -p = - d x E : | 
RN Í (24 a (2 J i ens = m-e) 


Cette relation montre que la perte de charge totale 4p, est proportionnelle à la 


— distance parcourue par le fluide. Dans ce cas, on dit que la perte de charge est 


réguliére. 


(T É 
Q, = d - | D° umm» C'est la formule de Poiseuille 
 128u u 2) 
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`> La perte de charge totale 4p, peut aussi être exprimée soit en fonction du débit 


> volumique, soit en fonction de la vitesse moyenne (débitante) de l'écoulement : 
D 
A — ( 
4 
4 


E TU Ар, ma 
=V 5 = р umm» 4 
Q =, 128u L da 
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> 4. Coefficient de pertes de charge 
. Généralement, la perte de charge est couramment exprimée en fonction de la 


m 


А 1 
pression cinétique (ou dynamique) de l'écoulement dans la conduite : Л pV 


_ Ainsi, la perte de charge totale s'exprime sous la forme : 


ap = EET | Жн: 2 J: e 


| р? D “ру? |2?" 
U 
64 u L 64u L L de 
pov. = praD D AD ^D 


Coefficient sans dimension 


. Donc, pour un écoulement /amínaire dans une conduite, on a : 


64 
À = CS I» est appelé coefficient de perte de charge réguliére (valable pour Re«2000) 
€ 
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> Dans le cas d’un écoulement dans une canalisation formée de plusieurs 
juxtapositions de conduites droites et de singularités, on peut alors généraliser 


l'équation de Bernoulli comme suite : 
A 


Singularités 


Singularités 


D, = PA* D 8 Za +S DVA 


: la t, 


L, 1 la d L 1 
w Ap =A, = оү? +A, 2 оү? +...+ A, =" — oV? +Pertes de charge singulières 
Р, “р, 5 P т 1 ^D, 5 P m2 " D, 5 P mn 


1 
(D, = ps +P 8 Ze DV 


> Of p, = p, + Ар, 


= 


C'est l'équation de Bernoulli généralisée 
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… Cette équation de Bernoulli généralisée reste valable quel que soit le régime de 


. l'écoulement (laminaire, transitoire ou turbulent). Par contre, l'expression de la 
perte de charge régulière A, doit être déterminé expérimentalement ou tiré 
* d'abaques (A, = 64/Re). 


. П пе reste alors qu'à étudier les pertes de charge singulières. 


> 5. Pertes de charge singulières 
llo considère la singularité suivante représentant un élargissement brusque, et 


à travers laquelle s'écoule un fluide incompressible. 


E: zones dites mortes ou 
zones de stagnation 


On considère la singularité suivante représentant un élargissement brusque, et à 


travers laquelle s'écoule un fluide incompressible. 


On choisit un volume de contrôle sur lequel on applique le théorème d'Euler. Pour 
cela, on supposera qu'on a: 
е - Pressions et vitesses uniformes sur chaque section (amont et aval), 
- Pertes de charge régulières négligeables sur le rétrécissement (sur L). 
Ainsi, la projection sur l'axe x, donne : 


Q, -(Vm 2 — Vm 1) = projection sur l'axe x, de la résultante des forces exercées sur le volume de contrôle : 


» poussée en amont : p,S, 

» contre - poussée en aval : - p,S; 

» poussée de la paroi verticale sur la partie stagnante du fluide : 
+ p(S; - S) où p = p, > (loi de l'hydrostatique) 


Comme le fluide est supposé incompressible, donc le débit massique se conserve: 


ny 


Comme le fluide est supposé incompressible, donc le débit massique se conserve: 
perte de 


1 20 Ï ee charge 
> Pit P Vm AS a 


p > Soit l'expression finale : 
] 2 
1 2 1 2 1 2 A J 1 
» +> оү? =p,+—pV2 + рү? | 1-2 | =p,+— рү? +— ру? 
| Pi 2 p 1 P 2 p m» 2 p m, | S, | P» 2 p m, 2 p m, 
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`» Ainsi, l'équation de Bernoulli généralisée s'écrit : 


Singularités 


Singularités 


Ce sont les coefficients de perte de charge singuliere qui sont 
associés à chaque singularité rencontrée au cours de l'écoulement. 


. Quelques singularités typiques : 
Elargissement brusque 


n S, 
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6. Pertes de charge en régime turbulent 
Lors d'un écoulement turbulent, les profils de vitesse dans une conduite 


cylindrique n'est plus parabolique. La turbulence de l'écoulement fait que les 
vitesses deviennent uniformisées sur une large partie d'une section droite 
donnée de la conduite. Ainsi, une brusque variation de vitesse au voisinage des 
parois est alors constatée. 


Ay Ar 


R 


Rë 


X1 
€ 0 


-R 
Écoulement turbulent 


Écoulement laminaire 
Re < 2000 umm À = 64/Re Re > 2000 шш) À = /(Re,c/D) 


L'équation de Bernoulli généralisée peut étre toujours utilisée, mais le coefficient 
de perte de charge 4 doit étre déterminé expérimentalement ou bien tiré 
d'abaques ou de lois empiriques. 
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Diagramme de MOODY 


Les tracés entirets 
0,080 ` ne Sont qu'indicatits ` 
) 


HH i 0,001 
= n racc s pi 0,0005 
SE went 


EARTH Lip —— Een EES 
LLL SER ъс SUED DRE 
Š 


ее 
Ф 
Ben 
2 
5 
si 
| `O 
[us 
Ф 
O) 
A 
G 
< 
. TO 
Ф 
© 
"|o 
L 
L 
Ф 
Q. 
|o 
© 
c 
„D 
o 
= 
Ф 
о 
o 


ESI 


7 0° Re 
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Détermination du diagramme de MOODY 


7. Détermination numérique du diagramme de Moody 
Étude de la perte de charge dans une conduite rugueuse 


Calcul du coefficient de perte de charge linéaire À dans une conduite 
On considère un écoulement d'eau dans une conduite de diamètre D. On veut 


=> calculer le coefficient de perte de charge linéaire À aux régimes laminaire et 


turbulent. 
En régime laminaire ( Re < 2000) la perte de charge est donnée par la formule de 
. Poiseuille : 
| 64 
E — 
Re (1) 


` ой Re est le nombre de Reynolds. 


En régime turbulent (4000 < Re < 105), ie coefficient de perte de charge linéaire 
À dépend à la fois du nombre de Reynolds Re et de la rugosité relative de la 
conduite z/D. 


> Pour une conduite lisse (£/D=0), cette perte de charge linéaire À est donnée par la 


formule de Blasius : 


À = —— (2) 
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Diagramme de MOODY 


Pour une conduite rugueuse (5/0 +0 ), la perte de charge linéaire À peut être 
calculée selon la relation (non linéaire) de Colebrook : 


E 2:31 
— + —— 
371.D Be Al? 


où À est le coefficient de perte de charge linéaire. 


A^ 2-2.log,, (3) 


1°) Écrire un algorithme (sous Matlab®, par exemple) pouvant évaluer à en régime 
laminaire et en régime turbulent (turbulent lisse et turbulent rugueux). 


Indications ` Pour l'estimé initial de Ay, on pourrait, par exemple, tracer la courbe 


- ҚА) en fonction de 2, ensuite repérer la valeur de À, pour accélérer la convergence. 


Si non, on pourrait toujours démarrer avec un estimé initial arbitraire et positif 
(valeur physique pour À). 
2°) Tracer sur le méme graphique (quadrillage) l'évolution de À en fonction de Re 


en échelle logarithmique (log-log). Pour l'écoulement laminaire, on choisit un pas 
de 20 pour Re et pour l'écoulement turbulent, on choisit un pas adaptatif pour Re 


(évolution de Re par palier). 
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Diagramme de MOODY 


3°) Comparer les résultats obtenus numériquement avec ceux données par le 
_ diagramme de Moody, étudié précédemment. 


On donne : 


e/D = de 10? à 10^ tous 107, 
de 2. 10^ à 10? tous 10^, 
de 2. 10? à 10? tous 10? 
de 2.10? à 7. 10? tous 10? 
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` Diagramme de MOODY 


Algorithme de calcul du coefficient de perte de charge régulière sous Matlab? 


L SKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


D + Pertes de charge régulière dans une conduite rugueuse : Ecoulements laminaire et Turbulent % 
gkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


b Close all;clear all;clc; 


Ë Régime laminaire 
D Rel = [0.1 1:1:10 20:10:100 200:100:2000]; 
11 = 64./Rel; 


> % Régime laminaire qui continue parfois jusqu'à Re=3500 
.Rel2 = [2010:10:3500]; 
L12 = 64./Rel2; 


figure (1); 
loglog(Rel,L1,Rel2,L12,':m','LineWidth',1.5);grid on; 
hold on; 


Ë Régime turbulent dans une conduite lisse : Loi de Blasius 
D Re2 = [4000:1000:1e8] ; 
L2 = 0.316*Re2.^(-0.25); 
loglog(Re2,12,'b','LineWidth',1.5); 
. &break 
% Régime turbulent dans une conduite rugueuse : Colebrook 
Re3 = [4000:1000:1е5 1.2e5:0.2e5:1e6 1.2е6:0.2е6:1е7 1.2e7:0.2e7:1e8]; 


% kl = [le-7 le-6 5e-6 1e-5:1e-5:1e-4 2e-4:1e-4:1e-3 1.2е-3:1е-3:1е-2 1.2e-2:1e-2:2e-1]; % kl = eps/D ===> Rugosité relative du tube; 
ki = [1e-5:1e-5:1e-4 2e-4:1e-4:1e-3 2e-3:1e-3:1e-2 2e-2:1e-2:7e-2]; % kl = eps/D ===> Rugosité relative du tube; 


————— 
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Diaqramme de MOODY 


4 ei 


Kä = 1/3.71; 
КЗ = 2.51./Re3; 


n = 20000; % Nombre maximal d'itérations 
Er = 1e-10; % Précision souhaitée 

% Relaxation 

w = 1.0; 


N1 = length(Kl); 
N2 = length (КЗ); 


| % Drapeau 
| s=0; 


> H£for i = 1:N1 
| for j = 1:N2 
Г clear x; 
x(1) = -10; % Valeur arbitraire initiale de Lamda 0 
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1 Diagramme de MOODY 


x for k = 1:n 

| f(k) = x(k) + 2*1ogl0(Kl(i).*k2 + k3(j)*x(k)); 

E dif (k) = 1 + (2/1og(10))*(k3(j)/(k1(i) *k2 + k3(j)*x(k))); 
x(k+1) = x(k) - f(k)/dif(k); 


x(k+1) = x(k) + w*(x(k+1) - x(k)); % relaxation 


if abs(x(k+1) - x(k)) <= Er 
x(k) = х(к+1); 


13 (j) = x(k)^(-2); % Changement de variable 
E N = k; s = 1; 
A break; 
4 end 
end 
end 


if (i == 1) | (i == 2) | (i == 11) | (i == 20) | (i == 29) | (i == 35) 
loglog(Re3,L3,'g','LineWidth',1.5); 
else 
loglog(Re3,L3,'r'); 
end 
end 
title('Evolution du coefficient de perte de charge en fonction du nombre de Reynolds'); 
axis([1e3 1e8 5e-3 0.1]); 
xlabel('Re'); ylabel('l','FontName','symbol'); 


ET 
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Diagramme de MOODY 


gtext('«--- Laminaire (Poiseuille) '); 


gtext('Rugosité relative'); 


gtext ('<--- Turbulent lisse (Blasius)'),gtext('e = 0','FontName','symbol'); 
gtext('e = 0.000001','FontName','symbol'); 

gtext('e = 0.00001','FontName','symbol'); 

gtext('e = 0.0001','FontName','symbol'); 

gtext('e = 0.001','FontName','symbol'); 

gtext('e = 0.01','FontName','symbol'); 

gtext('e = 0.07','FontName','symbol'); 

break; 

if s == 


figure (2) ;¿clear i; 
i=1:H1; 
plot(i,x.^-0.5) ;grid on; 


xlabel('Itérations');ylabel('Solutions Lamda'); 
title('Evolution de la solution : coefficient de perte de charge'); 
end 
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> 8. Puissance hydraulique 
La puissance hydraulique P, nécessaire pour véhiculer un fluide réel (visqueux) 
» dans une conduite comportant des singularités, est reliée au débit volumique Q, 


et à la chute de pression Ap par l'expression : |P, = др О, | 


Singularités 


Singularités 


1 1 L; 
Ap=pa-Pp = р 8 (zp -24)* 3P (v2, e À E EEN V: +D K; ei 
CN i i j 
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8. Exercice d'application 
On considére la conduite horizontale d'un pipe line de 150 mm de diametre 


transportant une huile en écoulement permanent, de masse volumique 

p = 900 kg/m? et de viscosité dynamique (и = 2x10? Pa.s, à la vitesse moyenne 
de 

V = 2 m/s). Le pipe line mesure 300 km de longueur et chaque station de 


pompage est distante de 30 km de la précédente. 


1) Quelle est la puissance hydraulique effective de chacune des stations de 


pompage ? 


2) Quelle est l'énergie totale dépensée pour transporter une tonne d'huile d'un 


bout à l'autre du pipe line ? 


Chapitre 5 


Analyse dimensionnelle et Similitude physique 
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>». 1. Introduction 
Quand un système étudié est trop complexe pour permettre une résolution 
1 complète des équations fondamentales, ou bien lorsque son comportement est 
chaotique, l'analyse dimensionnelle donne accès de facon simple à des relations 
entre les différentes grandeurs caractérisant ce systéme. 
= Ainsi, l'analyse dimensionnelle est une méthode qui permet d'obtenir des 
1 informations sur un phénomène à partir de l'hypothèse que ce phénomène puisse 
être décrit par une équation dimensionnellement correcte et liant certaines 
variables. 
. Avec un peu plus d'efforts, on obtient une solution partielle de presque tous les 
] problèmes, et non pas une solution complète. En effet, le mécanisme intime du 
| phénomène physique n'est pas révélé par cette analyse simple. 

Sur le plan pratique, l'analyse dimensionnelle d'un probléme permet : 

a) de réduire le nombre de variables et par suite facilite beaucoup l'exploitation 
des résultats, 

b) d'étudier sur une réplique à petite échelle (maquette) les performances d'une 
machine (prototype) à construire. Le passage, de résultats obtenus par 
l'analyse dimensionnelle sur la maquette, aux performances réelles du 
prototype, s'effectue à l'aide des regles de la similitude qui seront développées 
dans la suite de ce cours. 
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2. Analyse dimensionnelle 
2.1. Entités de la physique (grandeurs, unités et mesures) 


Tout raisonnement scientifique est basé sur le concept d'entités qui sont 
caractérisées par des effets observables. La force, la masse, le temps, 
l'accélération, la vitesse, la température, la charge, le courant électrique, etc., sont 
des entités observables. Chacune de ces entités possède une propriété appelée 
« grandeur ». Les « grandeurs » peuvent étre considérées comme des concepts 
abstraits. 
Par des procédés physiques, on peut établir une correspondance biunivoque 
entre les « grandeurs » d'une entité et une suite de nombres réels. On définit ainsi 
un systéme de mesure pour l'entité considérée. 
Pour mesurer la grandeur < G; > d'une entité (G;), le moyen le plus simple consiste 
à le comparer à une grandeur de méme espèce prise comme < unité >. 
La mesure s'exprime par un nombre réel « n; », qui indique combien de fois l'unité 
« u; > est comprise dans < G, >. L'écriture symbolique est donc : 

G; = n.u; (1) 
Si la grandeur Gi mesurée à l'aide de l'unité u;, s'exprime par un nombre (mesure 
п), la méme grandeur mesurée à l'aide d'une unité de mesure u’, s'exprimera par 
un nombre (mesure и”), tel que : 

П.и; = n'.u' (2) 
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2.2. Grandeurs et unités primaires et secondaires 
Suivant une convention trés largement respectée, trois entités peuvent étre 


considérées comme entités primaires ou principales. Ce sont : la longueur, la 
. masse et le temps (L,M,T). Ce choix ainsi que le nombre de ces entités, seront 
discutés ultérieurement. 


On attribue des unités arbitraires dites < unités primaires > à ces entités (u,, Um 
~ Uy). Par définition, les unités primaires n'ont aucune correspondance entre elles. 


Toutes les grandeurs physiques G, qui ont des caractéres définissables à partir 
. des grandeurs primaires par des relations monómes du type (3) sont des 
grandeurs secondaires ou dérivées. 


С, =F: М?Т” (3) 


| On peut aussi définir des unités dérivées permettant la mesure des grandeurs G, 


Soit, par exemple : i 


a l 
и, =иийи? m 
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II existe entre les mesures exprimées par des nombres réels n; des relations du 


type : 

п, = nn nZ (5) 
Les grandeurs "physiques Gi, dérivées des grandeurs primaires L, M et T sont 
symbolisées par des équations dites « aux dimensions ». 


(07) 1 
G, |=|" m”T”| g 
~ Nous allons limiter le nombre de grandeurs primaires à trois : longueur, masse et 
temps. Mais, on trouve aussi souvent : 
- En thermodynamique : la température беї la quantité de chaleur Q 
- En électricité : la charge électrique q. 
Dans ces deux cas, on peut montrer que si le choix de ces grandeurs comme 
grandeurs primaires peut s'avérer parfois utile, il n'est en aucun cas nécessaire. 


En effet, le premier principe de la thermodynamique dit qu'il y a équivalence entre 

~ énergie mécanique W et quantité de chaleur Q reçue ou cédée par le système en 
évolution cyclique, à condition de considérer Q comme un nombre pur. Dans ce 
cas, on peut écrire que : 


lo, |= р2.мт?| (7) 
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Quelques dimensions des unités usuelles en mécanique 


Entités — |Systémes (LM. 


Longueur 

Surface 

Volume 

Masse 

Temps 

Force 

Angle 

Vitesse angulaire 
Accélération angulaire 
Vitesse 

Accélération 

Masse volumique 
Poids 

Poids volumique 
Pression ou contrainte 
Moment d'une force 
Travail 

Puissance 

Quantité de mouvement 
Viscosité dynamique 
Viscosité cinématique 
Tension superficielle 
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2.3. Principe d 'homogénéité 
Les relations physiques que l'on traduit par des équations, font intervenir des 


opérations élémentaires de type somme, produit ou fonctions plus complexes qui 


| sont elles mêmes définies si et seulement si, la valeur donnée par la variable est 


un nombre pur. 
Ces relations doivent rester valables quelles que soient les unités adoptées en 
. restant toutefois dans un systéme cohérent. 
La découverte de la loi physique consiste, justement, à chercher la relation 
mathématique entre ~ G; : 
E J(G,,G;,Gs,G;.... ,G„) = (8) 
L'équation possédant la sasa kem d'étre indépendante des unités 
choisies pour mesurer les pice grandeurs, permet d'écrire : 
l KN. ;Nz,Nz,Ng,... DE (9) 
Remarque : 
Pour qu'une équation ait un sens de point de vue physique, l'homogénéité 
dimensionnelle est une condition nécessaire mais pas suffisante. 
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2.4. Grandeurs fondamentales 


Les systémes d'unité couramment utilisés reposent sur des unités primaires 
(L,M,T). Mais, on ne peut concevoir un système d'unités cohérent basé sur 


. d’autres grandeurs, dites fondamentales. 


Le nouveau système d'unités sera cohérent si une grandeur G; possède dans ce 
systéme fondamental, une unité de mesure définie sans ambigüité. 
. Soit, par exemple F, F, et F, les grandeurs choisies comme grandeurs 
` fondamentales. La grandeur G; doit avoir une dimension unique par rapport à F,, 
F, et F, et par suite, il doit être possible d'écrire. 


— Хд X;2 X;3 
G E |F F; e | (10) 
OU хү, X2; et X4; sont les coefficients inconnus recherchés. 


Dans un systéme primaire (L,M,T), les équations aux dimensions des grandeurs 
Gi, F,, F. et F, sont supposées connues : 


с, |= Ir Me T” ud 
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— C'est-a-dire : 


Ir, |- рмет“ | 
|F, |- IL M T| a2 
IF, |- |I M" T^" | 


. La substitution dans l'équation (10) de F1, Р, et F, (équation (12)) et l'identification 


| des exposants de L, M, T ainsi obtenus avec ceux de l'équation (11), conduit à la 
> résolution du système d'équations linéaires suivant : 


A X< T Xi FAAN = Ou 
Qt X; EE X; FAA = 055 (13) 
Az Xj T 0535 X; + 034 X4 = Or: 


Ce systeme doit admettre une solution unique non nulle, traduisant 


Q l'indépendance dimensionnelle des grandeurs F,, F, et Fz. 
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> La condition s'exprime par le déterminant non nul du système linéaire : 
Qj; Qi Ou 
А= 0, 05 0*0 (14) 
Gu @3 Gan 
De facon plus générale, soit p le rang de la matrice dimensionnelle du tableau des 
exposants dans le système primaire des n grandeurs intervenant dans le 


phénomene étudié (le rang d'une matrice est l'ordre le plus élevé d'un déterminant 
non nul formé dans cette matrice). 


II est donc possible de définir (n-p) produits sans dimensions, appelé souvent 


« groupements Æ >», constitués au moyen des p grandeurs fondamentales (ces 
grandeurs fondamentales sont celles dont les exposants dans le système 
primaire, forment un déterminant non nul). 


Soient : F, = G,, F, EG; ; Рз = Gs ; ... ; F, EG, 
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Les groupements z; associés aux grandeurs, dites dérivées, sont obtenues de la 


facon suivante : 
1 

Ch EE Ев Е | 

pod à 6e, 
Les exposants x, (p+1 < i < n et 1 < j < p) sont obtenus par la résolution du 
systeme d'équations (13) obtenu précédemment. 
On obtient alors (n - p) groupements sans dimensions. 
Vashy en 1890, puis Buckingham peu de temps aprés, avaient énoncé un 
théorème connu sous le nom de Vashy-Buckingham ou encore théorème z. 


(15) 


i 


`> Ce théorème dit, qu'une équation entre n grandeurs (éq. (8)) peut être mise sous 


la forme d'une équation entre (n-p) groupements sans dimensions, p étant le 
nombre maximum de grandeurs dimensionnellement indépendantes qui 
interviennent dans l'équation considérée. 
Ainsi, par application de l'équation (8), et en introduisant les groupements z, cette 
équation (8) devient grâce au théorème de Vashy-Buckingham : 

(тн, Tp+23 Tp43:---; Tn) = 0 (16) 
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. Remarque: 
> L'équation (16) peut aussi s'écrire de facon à particulariser un groupement т qui 
apparait alors comme une fonction des (n - (p-1)) autres : 


Tp+1 = V(n,.»; T3» `... Tn) (1 7) 


œ Dans tous les cas, l'analyse dimensionnelle ne peut donner aucune indication sur 
> la fonction F ou sur la fonction w. Seule, l'expérience peut permettre la 
formulation la plus complète d'une loi physique. 


4 
4 ^ 
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2.5. Comment poser un probléme en analyse dimensionnelle ? 
Lors de la phase préparatoire à l'étude d'un probléme par l'analyse 


dimensionnelle, il convient de : 
a) Récapituler les grandeurs actives, 


b) Choisir les grandeurs fondamentales, 
c) Former les groupements z suivant la technique exposée précédemment, 
d) Formuler la loi sous la forme : 


Amp Тр+2, 1рьз›· · Tn) = 0 (18) 
ou bien : 
VH = v(n,.. Тр+3›‹ een Тү» a... Tn) (19) 


e) Recourir à expérimentation pour préciser la loi physique. 


`> La loi, une fois formulée sous la forme (18) ou (19), n'est pas entierement définie. 


Car, l'analyse dimensionnelle ne peut pas faire apparaitre les constantes pures, 
les exposants, etc. Il y a lieu alors de recourir à l'expérience. 

On fait varier les groupements x les uns après les autres afin de déterminer leur 
degré d'intervention. 

L'expérience permet de mettre en évidence les grandeurs actives et/ou de faire 
apparaitre les grandeurs qui n'interviennent pas dans le phénomene étudié. 
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`» 2.7. Exemples 
Etude de la perte de pression motrice dans une conduite rectiligne 
Considérons un fluide de mase volumique р, de viscosité dynamique L, 


—  s'écoulant à la vitesse moyenne V, dans une conduite de diamètre D et de 


rugosité de surface e (hauteur statistique moyenne des rugosités). On souhaite 


trouver la loi de perte de pression motrice Ap sur une longueur AL de la 


= conduite. 
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| Solution 
Les grandeurs actives sont : D, e, AL, p, p, Ap et V. 
La matrice dimensionnelle est de rang 3. Elle s'écrit ainsi : 


L 


Prenons comme grandeurs fondamentales : D, p et V. 


= On vérifie dans un premier temps que ces trois grandeurs sont 


dimensionnellement indépendantes, puisque : 
1 -3 Ï 
А=0 1 0|=-1=0 
0 0 -I 
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On calcule ensuite les groupements x associés aux grandeurs : e, AL, u et Ap : 


Calculons par exemple le groupement Ap : 


. D’après le cours, on sait que : 
Л ap n D" p? V 
Les exposants x,, X, et x, sont obtenus par la résolution du système d'équations : 
x. =Ü 


(D.x +(—3).x, +(1).x, 2 -1 3 +24, +7; = 1 


(0).x, +(1).x, +(0).x; =1 
(0).x, +(0).x, + (-D.x, = -2 


| D'oü le groupement : 
) ЛАБ — — Soit : 
| D .p V 


=> 4x, =] On trouve alors : < x, = 1 
x, = 2 


=Y; ==? 
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. En procédant de la méme manière pour e, AL et p, on trouve : 


` Le théorème de Vashy-Buckingham s'écrit alors : 
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… L'expérience, ou plus simplement la connaissance élémentaire du phénomène 


étudié, permet de prévoir l'intervention du groupement AL/D comme facteur 


multiplicatif à la puissance 1. C'est-à-dire que la loi recherchée peut prendre la 


. . forme simplifiée suivante : 


À ce stade, l'analyse du probléme ne peut plus étre poursuivie, et il faut recourir à 


. lexpérience pour aller plus en avant dans la formulation de la loi. 
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> Remarque : 


Le groupement e/D s'appelle rugosité relative (s), et le groupement u/(D.o.V) est 


l'inverse du nombre de Reynolds Re. Il s'en suit que la loi peut s'écrire : 


AD AL 
- = ye, Re) 
PV D 
- Si Re < 2000 > l'écoulement est laminaire, 
— - Si Re > 4000 à 5000 > l'écoulement est turbulent, 
- Si 2000 « Re « 4000 à 5000 => l'écoulement est transitoire. 


La théorie des écoulements laminaires et les expériences effectuées en 
écoulement turbulent, confirment la forme de la loi que l'on écrit souvent : 


Ap AL 
ond 
PV D 
où À est le coefficient de perte de charge de la conduite, avec : 
64 
E Re > en écoulement laminaire, 


А= f(s,Re) > en écoulement turbulent. 
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Écoulement autour d'un objet 


On appelle force de trainée, la force 
exercée par un écoulement sur un objet, 


— dans la direction parallele à l'écoulement. 


CC š ou їч 


|> Nous allons étudier le cas d'une plaque - 
plane rectangulaire. Z 
1- Les variables du probléme sont : =s 


T, h, L, V, uet p _>n = 6 


T : force de traînée 

H : hauteur de la plaque 

L : largeur de la plaque 

V: vitesse de l'écoulement 
р: masse volumique du fluide 


aem 
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` 2- Équations aux dimensions : 


E " мг M "m 


FU us 
m së L WS A md. 
[]=[M'i'r'l=\M.lr3] pi 
[о]; IM pre |м] 


. 3- Nombre de produits x sans dimension : n- p = 6- 3 = 3 


_4- Choix de p = 3 variables primaires dimensionnellement indépendantes : 
> Par exemple h, pet V 


(Remarque : On ne peut pas choisir à la fois het L) 
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~ 5- Formulation des 3 produits x : ту = T/(he!.pP!. V1) 
Par combinaison des variables по = L/(h°2.p82.V"?) 
primaires et non primaire m, = u/(he3.p83, V»3) 
T 
Лу = ———— = Sans dimension 
h^ p” y^ 


ү ‚= " MIT Е (г MT? с“ fc: m'y’ (z T° |” » pia +3.8 -y, | y (1+0-6, +0) 72-0404 ) 


Ce qui se traduit par le système d'équations linéaires suivant : 


0=1-a, +3.) -7 - a +3.ß, - y, =-1 a, =2 
0=1+0- +0 ou: <f =1 ou encore : < f =1 


0=-2+0+0+7, 


=з 


dn = 


Soit : 
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—» Sans dimension 


kr ph y^ 
В z, PM = (E мот) ra^ (p ray intend atn ooon) 
. Ce qui se traduit par le systéme d'équations linéaires suivant : 


| 
) 0=l=0 +36 = у — E E a, =1 
0=0+0-, +0 ou : p; =0 ou encore : p; =0 


0=0+0+0+7, У =0 у, =Ü 


ч ж =É 


4 r 
d y 
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wu 
h^ p у? 


—» Sans dimension 


Ё =D.M° T° = (^ M'T! yr fc: Al" (z T” |” = Цаа) ule n р зана) 


Ce qui se traduit par le systéme d'équations linéaires suivant : 


4 pb 0--1-o,-43.f, - y, el 43.5, FY; = a, =1 
0=1+0- 2, +0 ou: 15, =1 ou encore : ; д, =1 
0=-1+0+0+7, y, =1 el 


À Soit : 
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[ 


. 6- Formulation des З produits x trouvés : 


T = f(h; LV,u, p) E = g[z,, m) 


| _ T L u 
oh h 


. avec: Z, 


Facteur de forme 


Nature de l'écoulement 
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> Illustration de l'intérêt de la méthode : 


Si T, est la trainée mesurée sur une plaque de dimension L, x h, quand elle est 
soumise à un écoulement de vitesse V;,, alors : 


T, = (E. 1 ) eu Re, -AV M 


2 
P М № h, Re, H, 
L'analyse dimensionnelle via le théorème de Vashy-Buckingham permet d'en 
. . déduire que pour une plaque de dimensions L, x h, telles que : 


m E dE 


gem de forme Va d'échelle hy deg 
еі done: g(L /h A/Re, )= g(L,/ h; 1/Ве,) etsi:p, =p, =p 
T, T, Vi 


- Ben o 1-5 > A 
oW) ovi Ж 
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— 3. Coefficients sans dimensions usuels 


Outre le nombre de Reynolds (pour définir le caractère laminaire ou turbulent d'un 
écoulement), il existe un certain nombre de grandeurs sans dimension qui 
. peuvent aussi caractériser la nature d'un écoulement : 


_ pV.L forces d'inertie 
Nombre de Reynolds : Ке = ——— 1 ЕССЕ СИ 


u 
L importance générale pour tout type d'écoulement. 


V forces d'inerti 
= Jre orces d'inertie 
Nombre de Froude : Bes: Ш> es s TC 


L importance pour les écoulements à surfaces libres. 


` Remarque : 


C'est la pesanteur qui est responsable de la forme de la surface libre : Plus Fr est 
grand, moins la surface libre a d'effet sur l'écoulement et inversement. 
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Nombre d' Euler : Eu = 


" mj. forces de pression 


forces d'inertie 


y? 
L importance s'il existe de grandes différences de 
pression au sein d'un écoulement. 


V forces d'inertie 
Nombre de Mach : Ma = — INZ > 


forces de compressibilité 


importance" pour les écoulements de fluides 
compressibles. 


p 1 
avec: a=./7.7.T = ly. = = —— 
P NEX 


et y= E Coefficient de compressibilité 
D \ CP /-isotherme. 


forces d'inertie locales 
Nombre de Strouhal : St = EA EE 


L forces d'inertie convectives 


importance pour les écoulements instationnaires. 
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Nombre de Weber ou : 


V : vitesse 
Nombre de Fourrier —— —, |і: longueur 
` p°: p : masse volumique 
C, : capacité calorifique 


Nombre de Nusselt À : conductivité thermique 


| t : temps 
Nombre de Sherwood K : coefficient de transfert 
de masse 
Nombre de Prandtl К D : Coefficient de diffusion 


u : viscosité dynamique 

a : coefficient de 
convection 

В : Coefficient de dilatation 

g 


Nombre de Peclet 


: accélération de 
pesanteur 
m: débit massique 
c : tension superficielle 


AT : variation de 


Nombre de Grashof Gr = Р température 


Nombre de Margoulis (ou nombre de Stanton 


Nombre de Graetz 


| 
] [Nombre de Prandt! ` ` 
|NombredePeclet ` 
| 
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> 4. Similitude dans les équations différentielles aux dérivées partielles 
Avant de procéder à l'analyse complète d'un écoulement, il convient tout d'abord 


de poser les hypotheses simplificatrices adéquates. 

. L'évaluation des différents coefficients sans dimension relatifs à l'écoulement 

| (nombre de Reynolds Re, nombre de Strouhal St, nombre de Froude FR, etc.) va 
permettre en effet de simplifier les équations à résoudre. 


.  Analysons ce que la composante verticale (selon l'axe z) de l'équation de Navier- 


Stokes, devienne aprés adimensionnalisation : 


Оуу Оуу Оуу Ор Ow Ow Ow 
ГВА Бота = VE FW. [= =H] =E = 4... на 02 O 

Ot дф 2 & Ox Oy д 

Introduisons le changement de variables suivant (sans dimension) : 

u*=u/V x*=x/L espi où L, V, p, et r sont des 
v*=v/V у*= y/L P =P! Po grandeurs caractéristiques 

* _ * _ t*=t/t d tè Studié 

w*-2wlV z*-zlL u systéme étudié. 
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—. Par extension, on obtient : 


9 190 
Ox L дх* 
9 1690 
Oy L Oy* 
ô 1060 
Oz L Oz* 


ü pV Ow* pV’ 
T ОР L 


0 1 © 
od L ôx? 
ó 1 © ô 1 ô 
o P ôy? Ot t Ot* 
0 ү © 
Oz? (DO dz 
Ainsi, la composante verticale de l'équation de Navier-Stokes devient : 

„ Ои“ Оуу * Оуу * 
Sp. +W*. 

ôx* ôy* Oz* 

__ Po Op* Jm [9w* 0w* дит De 
L @z* Ox? Gei Go" | 
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V Ow* PV Ow* =. ôw*  „ Ow* 
— + — > 


. 4 u * 


b т Or? L| Gr? i ôy* ôz* 


Op* иу [d'w* Ow* Owi 


. . C. x + C. xXx?) + x Nt 
L ô&z* L (O^ ду? Ge" 
2 
En divisant l'expression ci-dessus par le facteur 
J Оуу * Оуу * Оуу * BS 
. Tue Re SH w* 
Ot * Ce? 


, on trouve : 


== * I: д?у 
RESCH 


„= ү 


= L Fr = 
_ 


a 
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y 


~ D'où la forme finale de l'équation aux dérivées partielles : 
* 
# OW 
p 


](O^w* Ow* Owi 


Op * 
= —-Eu.—— + —. gU cur ER Ni 
Oz* Re! œ Oy Oz Fr 


D'où la forme finale de l'équation aux dérivées partielles : 


> Si Stest très faible Ə on peut négliger la dérivée instantanée et dans ce 
cas l'écoulement peut étre considéré comme stationnaire. 


> Si Eu est tres faible Ə on peut négliger le gradient de pression. 


> Si Reest très grand Ə on peut négliger la viscosité du fluide et l'assimiler 
à un fluide parfait (cas des très grandes vitesses d'écoulement). 


> SiFrest trés grand Ə on peut négliger les effets de la pesanteur. 
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A Similitude dans les équations différentielles 


Gráce au théoréme de Vashy-Buckingham, l'analyse dimensionnelle permet de 
résumer le comportement d'un système à une relation entre un nombre restreint 
… de grandeurs sans dimensions. 


Л, = ylz, л, Ty T) 


. Pour des systèmes complexes, la détermination de la fonction w n'est pas 
accessible que par des mesures expérimentales. Ainsi, lors de la mise au point 
d'un prototype, il est économiquement et pratiquement pertinent de procéder à 

_ ces mesures sur un modèle réduit : c'est la maquette. 


II faut pouvoir transposer les résultats obtenus sur la maquette à ceux que l'on 
obtiendra sur le prototype réel. 


Ainsi, on s'arrange pour respecter le maximum de similitudes entre la maquette et 
le prototype. 


E y 
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ь À titre d'illustration, nous reprenons l'exemple de la force de traînée T exercée par 
un écoulement sur une plaque plane : 


Maquette Prototype 
Т zu Hn | dy E H, | 
2 "e S бы 4 = Tr C 
ohh h, ok bh, Pp: VA h, h, p„V,.h, 
| | L 
. On commence tout d'abord par se fixer un facteur d'échelle : = -q-1/25 
Б | 
L, L 
On respecte ensuite le facteur de forme : P» E Bs 
m p 
` Si on utilise le même fluide, on а: P„=P, € H„=H, 
^» Respecter la similitude de Reynolds revient à écrire : 
V А VÀ h V 
Re, = Re es Pu nia Pt plo у V bh, =V, h, => V, =V = 
H, H, Р , а 
Т ү h 
Ainsi = 2 = (ECKER 
KE P V; AS Í 2 b 
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| Ce résultat, spécifique au probléme étudié, montre que par un choix approprié de 
Pla vitesse, la mesure expérimentale de Іа traînée sur une maquette 
géométriquement semblable donne directement la traínée à laquelle on doit 
| К s'attendre sur le prototype. 
A » En pratique, tout n'est pas simple, dans la mesure où l'on ne peut pas souvent 
. respecter simultanément toutes les similitudes. 
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` Comment garantir que les résultats obtenus sur une maquette soient 
… représentatifs des objets réels ? 
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Exemples de modeles réduits 


zi. 


SN 


inon (EDF) Érosion d'un barrage (CNR) 


t 


Maquette d'un sous marin (ONERA) 


а 


tte de ГАЗ80 (ONERA) 


|. Maque 


Écoulements dans les canaux découverts 


4 | 


Chapitre VII 198 


ÉCOULEMENTS DANS 
LES CANAUX DÉCOUVERTS 


> 1. Introduction 
|, M s'agit ici des écoulements pour lesquels il existe une surface libre 


. soumise généralement à la pression atmosphérique. Les canaux en 
| question peuvent être constitués aussi par les fleuves, rivières et canaux 

artificiels (par conséquent à ciel ouvert), que par certaines conduites 
`> fermées (tunnels, égouts, etc.) parcourues par des liquides n'occupant 


pas toute la section. 


rivière Canal artificiel 


199 


ÉCOULEMENTS DANS 
LES CANAUX DÉCOUVERTS 


L'écoulement dans les canaux est généralement produit par effet de la 
gravité, mais le mécanisme est plus compliqué que dans le cas des 
conduites en charge. Car la forme de la surface libre qui limite le 


domaine occupé par le fluide (liquide) est généralement inconnue. 


Dans la pratique, les écoulements sont rarement laminaires. La 


— turbulence se manifeste dés que le nombre de Reynolds dépasse la 


valeur de 2000 à 3000. Par ailleurs, le fluide qui s'écoule est 
habituellement de l'eau, et la plupart des mesures expérimentales 


existantes ont été effectuées avec ce liquide. 
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>» 2. Ecoulement permanent en moyenne 
2.1. Ecoulements uniformes 


2.1.1. Définition 
. L'écoulement dans un canal est considéré comme uniforme quand 


l'épaisseur du liquide le long de l'écoulement reste constante. Le plan de 
la surface libre est un plan incliné (d'angle o) parallèle aux génératrices, 
> comme le montre la figure ci-dessous. 

о 


Quand on coupe le canal par un plan perpendiculaire aux génératrices, la 
section correspondant au volume occupé par le liquide est appelée 
`> section mouillée du canal S. Le périmètre qui suit les parois solides du 
canal en contact avec le liquide est appelé périmètre mouillé Р. 


201 


ÉCOULEMENTS DANS 
LES CANAUX DÉCOUVERTS 


S 


. On définit le rayon hydraulique du canal par le rapport : Rg 
P 


r 


. Le diamètre hydraulique est donné par : Dy = 4.R„ 


Si on désigne par Q, le débit volumique du fluide à travers le canal, la 


> vitesse débitante U s'écrit : 


Q, 
S 


U — 


~ Le nombre de Reynolds Re qui sert à caractériser l'écoulement à surface 


libre à travers le canal est : 


po 4U Ry _ UDa 


V V 
Si a est l'angle formé par les génératrices du canal avec l'horizontale, la 


1 pente de ce dernier est définie alors par a-tg(o). 
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2.1.2. Cas des canaux en charge 


Si le canal est infiniment large, la vitesse est partout paralléle aux 
. génératrices de celui-ci. 


En première approximation, nous admettrons que la vitesse longitudinale 
| est constante et égale à U. De méme, nous pouvons admettre aussi que 
— les lignes de courant sont sensiblement rectilignes et paralleles aux 

génératrices du canal. Ainsi, la répartition de la pression est alors 

hydrostatique dans une section droite du canal. 
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— Le plan de référence étant le plan horizontal H passant par le point E. En 
un point M quelconque, on a : 


p+ p.g.h = Cte = p, 


~ En prenant la pression atmosphérique р, comme origine des pressions, 


et en comptant z positivement vers le bas (h = -z.cos(o)), on а: 

р 9.g.z.cos( a ) 
. Pour les canaux usuels, l'angle a est toujours trés faible, et par 
conséquent cos(a) = 1. D'oü : p = p.g.z. Cette expression est analogue а 


1 celle que l'on trouve pour un fluide au repos. 
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биг un élément de surface dS de la paroi solide du canal en contact avec 
le liquide, le fluide exerce la force de frottement 7,.dS qui varie avec la 


. position de l'élément considéré. Si 7,,, est la valeur moyenne de т, sur la 


surface mouillée, l'ensemble des forces de frottement s'exercant sur une 
longueur L du canal est égal à z7,„. P.L. 


— Les quantités 7,,, et 7; sont fonctions des conditions d'écoulement et de 


la nature de paroi (rugosité). Et comme pour le cas des conduites, on 

peut poser : 
~ U? 
Tom 7 Cy. p 

Les variations de C, sont analogues à celles dans les conduites. C, n'est fonction 

que du nombre de Reynolds Re et de la viscosité du fluide и (pour les régimes 


turbulents). 
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y 


~ Pour un écoulement permanent et uniforme, la relation liant les forces de 
> frottement aux parois sur une longueur L et les forces de gravité s'écrit : 
2 


U 
Tom P, L= C ,.p.— P. -L= p.g.S.L.sin( a ) 
> Comme l'angle a est trés faible, on peut alors remplacer sin(a) par tg(a)=i 


J. oü i est est la pente du canal, et on obtient la relation suivante : 


> Cette expression est connue sous le nom de formule de Chézy, où C est 


. le coefficient de Chézy (dimensions L'7.T") défini par : 


2. 
(=. 72 
C, 
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2.1.3. Calcul du coefficient de Chezy : Formule de Bazin 


Le coefficient de frottement peut étre exprimé en fonction du coefficient 
. de perte de charge linéaire 2 (dépendant du nombre de Reynolds Re еї de 


— la rugosité des parois lorsque Re > 3000) раг: 


б 


f 


> D'oü le coefficient de Chézy, devient : 


Cette relation est assez bien vérifiée par l'expérience tant que la largeur 
du canal n'est pas trop grande par rapport à sa profondeur. 


207 


ÉCOULEMENTS DANS 
LES CANAUX DÉCOUVERTS 


Lestracés entirets 
0080 ` ne sont qu'indicatits ` 
А 


0050 
oos E X | 
0,040 E SS 
0,035 


TT 
RSR 


| 


Eu] 
LL 


0,020 2 
wi: __ Hi 
due H 


Se | ke ЭҢ E3 SE BS LE Е AS Л БЕ 


|| LU 

EE FIRE 
EE 

= 


LE 0,0001 
—=== 


F FIM 4 
7 


10 10 


ее 
pum 
“Ф 
Е 
О 
`D 
x 
D 
љ 
G 
= 
(S) 
"© 
Ф 
C 
љ 
Ф 
Q. 
"© 
Amt 
C 
E 
o 
s= 
Ф 
° 
O 


HE 


208 


ÉCOULEMENTS DANS 
LES CANAUX DÉCOUVERTS 


= D'autres formules, d’origine empirique, ont été proposées pour le calcul 


de C. Parmi ces dernières, la formule de Bazin (ci-dessous) est 
. couramment utilisée (surtout en France), où le coefficient C est exprimé 
en m/s. 


. Le coefficient y prend les valeurs suivantes : 

+ у= 0,06 pour des parois très unies : ciment lissé, bois raboté 
* y = 0,16 pour des parois trés unies : planche, brique, pierres de taille 
* y = 0,46 pour des parois en maçonnerie de moellons 

* y = 0,85 pour des parois en terre régulière 

* y z 1,30 pour des parois en terre ordinaire 

* у= 1,75 pour des parois herbées, des fonds de galets. 
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2.1.4. Débitance d'un canal 
La détermination du débit d'un canal de forme 
. géométrique quelconque et de nature de parois 
| donnée, est un problème important dans la 
pratique. Elle s’effectue aisément si on connait la 
. section droite mouilée du canal S et la vitesse 
débitante U. 


Mais comme la section S et le périmètre mouillé P, dépendent de la 

. profondeur у du liquide dans le canal, le coefficient de Chézy C, qui ne 
dépend que de H,, va dépendre aussi uniquement de y. Ainsi, le débit 
volumique Q,, sécrit alors sous la forme : 


OÙ K = S.CAR, est une fonction de y appelée débitance du canal. 


H 
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> Si on examine l'évolution de K(y) en fonction de 


de K(y). 


la hauteur y, on obtient la courbe ci-joint. La 
quantité О, / Vi est tout simplement l'expression 


K(y) 


Les vitesses qu'on peut admettre dans un canal ne peuvent étre 


> quelconques à cause des risques d'érosion des parois. Elle atteignent 


0,20 m/s pour les terres argileuses, 1 m/s pour les graviers et 4 à 5 m/s 
pour les roches dures. 


. Les cours d'eau natuels ont généralement un lit assez inégal qui rend 


parfois assez difficile l'application des formules établies précédemment. 
| Cependant, on peut souvent assimiler leur section à un rectangle de 
1 largeur /, assez grande devant la profondeur y. Dans ce cas, 5-/.у, et 
5/Р,=у./1+2у) # у, et on obient: p _ [У 
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2.2. Ecoulements graduellement variés 


. 2.2.1. Définition "UU 
On dit qu'un ecouiement est graduelement varié si la pente de la surface 


| libre et celle du canal, ainsi que la forme de la section transversale, 


varient suffisamment lentement pour qu'on puisse admettre que 
l'écoulement ait un caractere régulier assez proche du régime uniforme le 


> long d'un petit trancon. En première approximation, on peut admettre que 


dans une méme section droite, la vitesse est constante, et égale à la 
vitesse débitante (UZQ/S), et que la répartition des pressions est 


`» hydrostatique. 


Lorsqu'on barre, par exemple, un canal en écoulement uniforme par un 
| barrage sulélevant le niveau en une section, l'effet s'en fait ressentir loin 
`` en amont et se manifeste par des variations de la profondeur d'eau et la 
pente de la surface libre. La variation de celle-ci par rapport à celle qu'on 


C avait en régime uniforme s'appelle le remous. 
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Lorsque la гуа аи canal est faible, la ligne piézométrique est 
pratiquement confondue avec la surface libre. La ligne de charge est 
. toujours descendante (dans le sens de l'écoulement). 

Si on se référe au canal et non ln 


| pas à un plan horizontal de — .— Ой ire de charge 
` référence (arbitraire), la cote de Uz _ eg — 
la ligne de charge par rapport T B E k i 
. au fond F du canal est : pé FES RR / C" Е 
| C2 PU 
I "Side; 
`> H, est appelée charge spécifique dans la Plan de référence ` 22 


section. C'est une quantité généralement différente de la charge H, mais 


Е elle peut servir à caractériser l'énergie totale du liquide. 
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2.2.3. Nombre de Froude 
Le nombre de Froude (de l'hydrodynamicien anglais William Froude), est 


un nombre adimensionnel qui caractérise dans un fluide l'importance 


`> relative des forces liées а la vitesse et à la force de pesanteur. Ce nombre 


apparait essentiellement dans les phénomènes à surface libre, en 
particulier dans les études de cours d'eau, de barrages, de ports et de 
= navires (architecture navale). 


En dynamique des fluides il fait partie des trois nombres sans dimension 


`> les plus utilisés (Re, M et Fr). Il caractérise l'importance de la pesanteur, 


alors que le nombre de Reynolds He prend en compte la viscosité et le 
nombre de Mach M la compressibilité. Ce nombre est défini par : 


où U est la vitesse de l'écoulement, g est l'accélération 
de la pesanteur (9,81 m.s) et h est la hauteur d'eau. 
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Pour un cours d'eau donné, un méme débit peut être obtenu de deux 
facons différentes : 
e Si Fr < 1 : Le régime de l'écoulement est dit “fluvial”, avec une forte 
hauteur d'eau et une faible vitesse (équivalent d'un écoulement 
subsonique). Ce régime est "piloté par l'aval" : le comportement des 
particules en mouvement est contraint par celles qui les précedent. 


e Si Fr > 1 : Le régime de l'écoulement est dit “torrentiel”, avec une 
faible hauteur d'eau et une forte vitesse (équivalent d'un régime 
supersonique). Dans ce régime, le fluide est "tiré" par les forces qui le 
meuvent (la gravité le plus souvent), sans que la masse de fluide en 
aval soit une gène. 
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Dans les deux solutions (Fr « 1 et Fr » 1), la hauteur d'eau et la vitesse 
sont déterminées suivant le nombre de Froude et le débit. Mais les 
solutions ne se calculent pas de la méme facon. La détermination du 
nombre de Froude est donc un préalable au calcul. 


La transition du régime torrentiel au régime fluvial provoque un ressaut 
hydraulique oü la hauteur d'eau s'accroit brusquement. 
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2.2.4. Ecoulement à débit Q, constant, y étant variable 
— La cote du fond étant donnée, il est plus simple d'étudier la charge 
comptée à partir du fond, c'est ce que nous l'avons appelée la charge 
| spécifique H. 
. Cherchons comment varie H, quand, dans un canal, s'écoule un débit 
Q,=S.U constant : , 
H =y +- #? 
m" 2.0.5 


ЕР 
= 


S est une fonction de y, qui varie dans le 

J même sens que y. Quand y est nulle, S est 

nulle et Hs devient infinie. Quand y est 

| infinie, S l'est aussi, mais H, l'est 

“également. Comme Hs possède des 

| valeurs finies, ses variations en fonction 
CN . de y doivent passer par un minimum. 
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— Ce minimum de H, est obtenu pour : 


ан, d( Q | , Œ ds_ 
` dy dy| 2.2.8" gS? dy 
Comme dS/dy = l (largeur de la surface libre), donc 
Q; 45 O/ 
gS? dy g.S? 
Soit, d’après la minimisation de dH./dy : 
2 
l 
0 
| g.5 
»La valeur 


de y correspondante est 
. appelée profondeur critique y. Elle se 
rapporte à la section critique S.. 
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lI Ainsi, à toute valeur suffisante de H., correspond donc deux valeurs de y, 
l'une en M (y > y.) pour laquelle le régime est fluvial, l'autre en N (y < y.) 


> pour laquelle le régime est torrentiel. 


En régime permanent, le méme débit Q, |4 
fluvial „B 


=~ peut donc s'écouler dans un canal de : 
torrentiel 


deux manières différentes, tout en ayant 
.. la méme énergie totale : 


* Régime fulvial : la vitesse est faible 
mais la profondeur de l'eau est grande, 
* Régime torrentiel : La vitesse est forte, 


mais la profondeur de l'eau est faible. 
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Relation entre dH./dy et le nombre de Froude Fr : 
1 Ona: 


2 2 
e Al + |+ =, йз _у Qd 


gs dy 
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> 3. Etude de la charge hydraulique en écoulement turbulent 


>La charge hydraulique est définie à partir des équations de Naviers- 
~ Stokes turbulentes. On moyenne ensuite cette charge sur la section d'un 


» écoulement à surface libre dans un canal. 


Ew 
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E 


harge hydraulique des écoulements à surface libre 


(5 
HES Y= Gs + Z (s) + h(s) + (в)? 


1 
А(5) 


H(s2) — H(s1) = — [ J(s) ds avec J(s) = 


fre 
A 
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A Ressauts hydrauliques 


| 
» 
d 


Le tracé de la charge spécifique permet de décrire l'écoulement 


… stationnaire sur un obstacle d'extension finie et de faible pente. Le tracé 


». de l'impulsion permet de décrire les ressauts stationnaires. 
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Équations de Navier-Stokes en moyenne de Reynolds : 


dvU=0 e 2+0. grad Hess Prado tu DU iem 
: R 


rot et F — —grad (g z) 


[Cx v AU div R) oi 
L Ot m 


il 
et {== pane 
pg 2g 


Integration de Mı à Mo 


H(M2) = НОМ) — f, (35g +). dM 


avec J= 1(— v AU + div R) 
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@ Section : A e Rayon hydraulique : Ry = A/P 


e Périmètre mouillé : P e Diamètre hydraulique : Dy = 4 Ry 


e Hauteur hydraulique : hy = A/ L 


e Largeur mirroir : L 


e On suppose souvent hy ~ h 
ә Section rectangulaire avec L> h: Ry ~ h 


e Section en demi-cercle de rayon R : Ry = R/2 
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Équations d'équilibre avec œa > 1 : 


d dH U? 
Sw р avec Q= UA, —— —J avec Н = д. +h+ Z += 
ds ds P E 28 


Pente / = —– 02; / ds : 


Si dA ~ Ldh et A/L = hj ~h: 


AUA) dU „ y dA dh _ Ç dh _ I-J 
ds `` de dh ds ` de 1— Fr? 


oü le nombre de Froude est 


dH _ dh , _ U dU _ 
ds ds gds Fr — Uf gh 
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I 
' 
L 


we e а mn l. mp mm I mme me em ФА ЫШ", 1 
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Hauteur critique he dans le cas L constante et A/L = hy ~ h : 


dh I — dZ; / ds 


ds 1—Fr? 1 (h/h) 


2 
c s LR i кок 
a = и, 


On a donc:  Z[q,h(s)] = Ho — 7[=) 


2 
H= NA OT SIR d => ГЕТ sma 
pg 2 8 P g 
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Solutions de Efq, h(s)| = Ho — Zr(s) 


€ Dän | | EE 
LL === 


Écoulement critique : Fr 


e Sous critique : 
Н => fi c= Fr< 1 


e Super critique : 
et => ЕР 3 1 


hı et h2 sont des hauteurs 


LU EE ul JE NE. conjugées pour E(q, h) 


EMENTS DANS 
ча - 11 VAI i ә 


| DÉCOUVERTS 
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Dur. — “a ef" Mal Se 
"we, 
X > > ien 


= 4 "rm V 
x “к= 


| - 
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y 


> 5. Courbes de remous 


Pour obtenir des courbes, on paramétrise le frottement, qui est 
~ proportionnel à la perte de charge, à l'aide d'un coefficient de Strickler, et 


| on trace les courbes de remous des écoulements stationnaires. 


La rugosité du canal intervient donc dans le calcul du coefficient de 


Chézy C, gráce à la loi de Manning-Strickler : 


C — KR" 


. oü K est appelé le coefficient de Manning-Strickler. 
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b Équations de Navier-Stokes turbulentes et perte de charge : 


> 1 1 
J = —(—v AU+divR)= —— div (т avec т. = p(2vd —R 
B el ao F. z,=p(2vd-R) 


J=} ff 2 e de 
A JJA 

SE dl 
fai i „ês mon 


IRAC EEN da 0 


H Relation entre la perte de charge linéique J et le frottement т, : 


F; =p E Ry J avec Рн = A/P 
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Paramétrisation du frottement т, = pg Ry J : 


U? 


т = = A(Re, Ru) p Ш? => J — A(Re, Ru) 5 
avec Re— U Dy/v et Ru = k;/Dy 


Ja let Ву = h => 


U = К, hš \/sin y 


235 


ÉCOULEMENTS DANS 
LES CANAUX DÉCOUVERTS 


Équation stationnaire pour h avec frottement J : 


dh _ | 11— (h/hs)- 19/3 
uu c sis (h/h.) 3 
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1 — (h/h,)- 19? 
1 (h/he) ^? 


а — 
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6. Exercice d'application : Etude d'un diversoir à seuil large 

Le seuil AB d'un déversoir de largeur L (perpendiculaire au plan de figure ci- 

dessous) est suffisamment long (AB) pour que la lame d'eau qui traverse la 
. section S puisse être considérée comme constituée de filets parallèles horizontaux 

s'écoulant à vitesse uniforme. La pression atmosphérique p,,, est également 

considérée comme uniforme. 


1) Calculer en fonction de la baisse de niveau y le débit qui traverse la section S. 
2) Pour quelle valeur de y le débit est-il maximum ? 
3) Exprimer ce débit maximum si la vitesse V, est négligeable. 


Chapitre 7 


Ecoulements rampants laminaires et lubrification hydrodynamique 
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ÉCOULEMENTS RAMPANTS LAMINAIRES 
ET LUBRIFICATION HYDRODYNAMIQUE 


1. Introduction 
Un fuide newtonien incompressible est dit en écoulement de Stokes ou 


en écoulement rampant, si l'on peut negliger les termes inertiels dans 
l'equation de Navier-Stokes devant le terme visqueux diffusif. Le terme 
non linéaire (dit aussi « convectif ») peut étre considéré comme étant 
| négligeable si le nombre de Reynolds He reste inférieur à 1. 

Dans се chapitre, on considére que les écoulements sont continus et 
permanents. Les fluides visqueux étant supposés isovolumes et non 


— pesants (dans le cas faibles epaisseurs de films liquides). On considère 


également que ces écoulements sont plans (bidimensionnels). 
De tels écoulements se présentent dans la /ubrification par liquide, du 


Э moins quand l'épaisseur du film liquide est suffisante pour que ce film 


puisse étre considéré comme un milieu continu. On dit alors qu'on a 


Q réalisé une < lubrification hydrodynamique >. 
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Nombre de Reynolds : Re 


Forces d'inertie du fluide 


Forces visqueuses 


OM Energie cinétique du fluide 


= Energie dissipée par frottement 


Temps de transport de QDM par diffusion 


— Temps de transport de ОРМ par convection 


B 
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Classification des écoulements : Régimes 


laminaire turbulent 


Re ««1 


Еа 


mi 


100 - 1000 He 
Filets fluides Mouvement 
paralléles désordonné 
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>» 2. Ecoulements entre deux plans parallèles 
Si l'écoulement est parallèle à une direction donnée ox, les lignes de 


courant sont des droites parallèles. Le calcul s'effectue comme pour le 


=> cas d'un tube cylindrique, et dans le cas d'un écoulement permanent, on 


trouve : У 
UP T — ufy : 
| h-2yoi — кз шлш = >> 
ой u ne dépend que de y et p ne -Yo bm 
> dépend que de х; ce qui conduit à écrire дие: 
dp — ( = 
dx 


> La pression varie linéairement le long de l'écoulement, et pour une 
_ conduite de longueur /, опа : 


— A 
GE p, ^P 


l l 
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— Ainsi, l'intégration de l'équation différentielle du second ordre : 
2 
du Іар -a 


d шах A 


` fournit la solution générale suivante : 


EE 
u 2 


où b et c sont des constantes d'intégration à déterminer par les 
conditions aux limites. 


>» 2.1. Ecoulements plans de Poiseuille 


Les deux parois sont au repos. Donc u = 0 à yzzy,. On trouve : 
2 


a >» y 
u( y= |е 
2и Уо ay 
a vitesse maximale est atteinte au centre (en y=0). Elle vaut :u( y ) = ; "i 
u 
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— La répartition des vitesses est donc parabolique et va dépendre du signe 
du gradient de pression : 


typo n 
` C \\ U d 
pay À Plans immobiles 3 = 
Уо РРР 
У. ү 
+Vo 4 
p-o, L. r=" _ _ _ Plan supérieur en w o 
=<Yo\ D U(y) mouvement de translation dx 
Yo Z 
У. ү 
+Vo à SEY ————— 
l ф===—— Ari d 
Aal Se Hl A Plansupéreuren ` 5 ©, 0 
ОГ —— mouvement de translation dx 
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Exercice d'application 


On considère l'écoulement d'un fluide de viscosité dynamique и entre 


> deux plans horizontaux distants de ^, l'un étant animé d'une vitesse V 


parallèlement à l'autre fixe, et soumis à un gradient de pression motrice 
^ 
longitudinal d p/ dx. 
1) Établir le profil de vitesse entre les deux plans, et étudier les cas suivants : 


Zp_o 
dx 
2727 
dx 
ар 0 
ах 


2) Donner l'expression du débit linéique q pour une largeur f. 
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a, 2.2. Ecoulements de Couette 


d | | | 
Dans се cas, on a P = —q = 0, une paroi est immobile, l'autre en mouve- 


ment avec la vitesse uniforme V. Le profil de vitesse devient alors : 


ууз E s. 
2 Yo 


La répartition des vitesses est donc linéaire. Elle correspond à l'exemple 


` que nous avions étudié dans l'exercice précédent. 


. Pour l'écoulement de Couette généralisé, le gradient de vitesse n'est pas 
toujours nul. Il peut être positif ou négatif. Et c'est aussi l'exemple qui a 


été étudié précédemment. 
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A. Ecoulements de Stokes ou écoulements rampants 
Lorsqu'un fluide visqueux s'écoule /entement en un lieu étroit ou autour 


d'un petit objet, les effets visqueux dominent sur les effets inertiels. Son 
écoulement est alors appelé écoulement de Stokes (ou écoulement 
rampant ; et on parle parfois de fluide de Stokes par opposition à fluide 


| parfait). || est en effet régi par une version simplifiée de l'équation de 


Navier-Stokes : C'est l'équation de Stokes, dans laquelle les termes 
inertiels sont absents. Le nombre de Reynolds Re mesure le poids relatif 


> des termes visqueux et inertiel dans l'équation de Navier-Stokes. 


L'écoulement de Stokes correspond ainsi à un faible nombre de Reynolds 
Re (Re << 1). 


`> L'équation de Stokes permet en particulier de décrire les écoulements de 


liquide dans les dispositifs microfluidiques. Les écoulements de Couette 
et de Poiseuille sont aussi décrits par cette équation. 
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| Effets visqueux 
| Diffusion de la quantité 
» de mouvement 


- ) H u V ^v 


O(uUJ L?) 


Effets d'inertie 
Transport de la quantité 
de mouvement 


4 z 
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Théorie de la lubrification hydrodynamique 


Hypothése géométrique de base + Valable dans la 
zone centrale uniquement en théorie 
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II —» Conservation la masse 


Qu "- Әр 4 ‚^ч 1 ‚^ч 
ou — +pu— = Lu — - 
дт e ду Or á ag? ` ' 


ду? 


m дь — Ov др = dv " ‚2° 
Ред, Ti Oy | dY Р Өг? "aya 


lee quantités de mouvement 


Que deviennent ces équations ? 


| h < L L = 
i è, D 
^ 
hol d 2 LR M 
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Ou du 


ti 


1 Yaw 
dy? 


L P 
Or 2 


Longueur horizontale caractéristique : L 
Longueur verticale caractérisitique : h 
Vitesse horizontale caractéristique : U 
| Comment choisir une 
vitesse verticale 


caractéristique ? 
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Quand peut-on négliger les 
termes d'inertie ? 


O(pVU/h) O(nU/ L*) < O(uU/h?) 


Hypothése de lubrification : 
Ecoulements rampants 


UL _ pUL M 
orces d'inertie | | pU*/L р y und 


mam e TIR — L? 
| Orces visqueuses uU /h H 
I "L 


Rer, 
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Et l'autre équation ? 


O(pU?h/ L?) 


O(nUh/L?) < O(uU/Lh) 


On obtient la 
méme condition... 


gut _ р M 


orces visqueuses uU / Lh p 1? 
MA 


Rer, 
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Et la pression ? 


p(r,y)— po = |р(т,0) — po[- E 
rd 
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Equations 
de Reynolds (1889) 


Film fluide mince 


Qu | Qv 


Or ду 


h < L 


0 = 


Hypothèse de lubrification : 
Ecoulements rampants 


pUL h? 
= < 1 


и І? 
SÉ 


Rer 
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Chapitre 8 


Aérodynamique externe en fluide idéal 


К f 
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AÉRODYNAMIQUE EXTERNE 
EN FLUIDE IDÉAL 


>» 1. Introduction 


Par définition, l'Aérodynamique est une partie de la physique qui s'intéresse aux 
q phénoménes accompagnant tout mouvement relatif entre un corps et l'air ou il 
` baigne. 


Les *Phénoménes" particulièrement intéressants pour l'ingénieur sont : 


+ Efforts subis par des véhicules (avions, automobiles, trains...) ou des 
projectiles (fusées, missiles, etc.) en mouvement dans l'atmosphère : 


> AERODYNAMIQUE EXTERNE 


+ Etude des systémes de propulsion de ces véhicules ou projectiles qui 
exploitent les propriétés de l'air ingéré au niveau des entrées d'air : 


> AERODYNAMIQUE INTERNE 
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Les ojectifs essentiels recherchés sont : 


Détermination de la trainée et de la portance 
occasionnées 
par l'écoulement externe autour des coprs, 


Comprendre de concept de la couche limite sur la surface 
de l'obstacle placé dans l'écoulement, 


Pouvoir calculer de la couche limite (laminaire ou 
turbulente) sur la surface de l'obstacle. 
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xÉcoulements externes : 


1- Portance et trainée 

2- Couche limite laminaire 

3- Transition laminaire/turbulent 
4- Couche limite turbulente 


5- Coefficients de trainée et de portance expérimentaux 
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~ Portance et traînée : 


Lorsqu'un corps se trouve complètement ou partiellement dans 
. un écoulement externe, il subit deux types de forces : 


=> + Des forces perpendiculaires NS 
à la surface du corps qui sont y mE | 


les forces de pression. 


T Distribution de 


p 
` — PA la contrainte de 
* Des forces tangentes à la SE — cisaillement 
surface du corps qui sont LA f VA. 


A. 
e А — ` — 
des forces de cisaillement. N — — — 
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Trainée et Portance : 


Les forces de pression et de cisaillement éveloppées sur la surface d'un 

` corps sont décomposées en deux composantes : 

1- Force parallèle à l'écoulement incident appelée force de trainée 
(généralement notée D, Drag). 

2- Force perpendiculaire à l'écoulement incident appelée force de 


portance (généralement notée L, Lift). 
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— Trainée et Portance : 


. Hypothèse 1 
«Ww Lorsqu'un corps se déplace dans fluide 


stationnaire (mouvement d'un avion dans 
une stratosphère stationnaire), ou d'un 


écoulement autour d'un corps fixe. 


=> On fixe un repère lié au corps et on 


considère que le fluide arrive avec une 


Corps 


vitesse incidente constante V. 


c» Écoulement incident uniforme et constant. 


1111111111 
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— Trainée et Portance : 
» Hypothése 2 


La structure de l'écoulement est classée en trois 


. catégories générales selon la forme du corps : 


. а) Bidimensionnel (2-0), lorsque le corps est | 
infiniment long avec une section longitudinale — — + 


parallèle à l'écoulement constante 


b) Axisymétrique : lorsque le corps a un axe de А — TENET 
AN 
symétrie dans la direction de l'écoulement z 


c) Tridimensionnel (3-D) 
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~ Trainée et Portance : 


Ces coefficients peuvent étre calculés à partir des expériences en se basant sur 


l'analyse dimensionnelle ou aussi par des méthodes numériques. 
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Classification des écoulements 


p? 


1*' Principe de la thermodynamique 


Lois de conservation : 


- de la masse, 
- de la quantité de mouvement 
- de l'énergie 


Ecoulement incompressible 


peche ERE нна зак кае вина кАй а. 


í [Ecoulement compressible| : 
E : Nombre de 


M=0 M = 0,3 : M=1 M>1 : MachM 
Régime subsonique : Régime transonique Régime : 
A supersonique 
Couche limite : 
: : Nombre de 
Re > 0 : :Ве 2 o» Reynolds Re 
: Régime laminaire Régime turbulent 


` 
SCTTETETETERTTTLTTEEITETLTLTTTTTTTTELTTTTTTTTTTTTTTTTTTTTTTTETTTTTTT TT T TTTTTTTTTTTTTIT LM 


АШ HY Ü A = 
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Positions schématiques de quelques domaines 
d'applications dans le plan Mach-Reynolds 


Nombre de 
eynolds Re _ 


Ве > œ ~ | | DEDE 
Véhicules de série Véhicules de course Avions de кипе MESES е 
e Chasse Fusées 
lignes 
Approche Croisiere 
@ 
`O 
© 
`O 
o 
o 
Been 
а. 
Compressible 
Re > 0 P Nombre de 


0 Incompressible 1 Mach M 


A A 
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Outils d'analyse des écoulements 


+ Intérét majeur Ə Prédiction Computational Fluids Dynamics (CFD) 


+ Principaux inconvénients > + Intérêt majeur > Souple d'utilisation et 
méthode tres cóuteuse, limitations | peu coüteuse(en général) 


Dépendance des résultats vis-à-vis de 
la modélisationet parfois de la méthode 
numérique > Pas toujours prédictive 


Full-span 
model s 


` Г, Г 
10 — Leg | Mx E... 
2 Cor«entional Wind Iuprele — 


Ü 02 04 06 08 10 12 


Données expérimentales (soufflerie Calculs Boeing Calculs ONERA 


2 


o 


8 


AÉRODYNAMIQUE EXTERNE 
EN FLUIDE IDEAL 


AM A 


Aérodynamique 
expérimentale 


Aérodynamique 
fondamentale 


Aérodynamique 
numérique 
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Modeles de Fluide 
Approche microscopique : 
On s'intéresse au comportement de chacune des molécules du fluide 
= Théorie Cinétique des Gaz 
Approche macroscopique : 
On définit un comportement moyen des molécules comprises dans un 
volume élémentaire (particule fluide) = Mécanique des Milieux Continus 


> - Approche pertinente selon la densité moléculaire (caractérisée par le libre 


parcours moyen Л): 
v Cinétique des gaz > Spatial 
v Milieu Continu > L'essentiel des problèmes en aéronautique, 
domaine automobile et des procédés, etc. 
- Paramètre distinctif : nombre de Knudsen = Л / L «0,01 el > MMC 
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» K Coefficients de pression 
_ Si p est la pression statique en un point de la paroi du corps, on définit ainsi le 
coefficient de pression locale par : 


~ C'est un coefficient sans dimensions, où Pœ, V. et po désignent toujours les 
.. conditions à l'infini. 

Pour un fluide incompressible et parfait, le théoréme de Bernoulli, appliqué entre 
=> un point situé à l'infini et un point situé à la frontière de la couche limite 


(conditions p, p. et V), donne l'expresion suivante : 


1 1 
+= p Vi =p+=p„V' 


Et le coefficient de pression locale devient alors ainsi : 
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`> Cette relation permet de définir les pressions, connaissant le champ des vitesses. 


Pour un fluide réel, lorsqu'un sillage se forme à l'aval du corps, la relation 
=> précédente, reste valable avec une trés bonne approximation, tant que l'on ne peut 

négliger la compressibilité jusqu'au point de décollement D. 

Mais dans le sillage, il ne peut étre question d'utiliser cette formule pour la 
`> détermination de la pression en un point de la paroi. 
D 


Boundary layer p 
separation 
location 


7 zu Леа: o ага сь 


Boundary layer 
D fluid has kinetic 
energy deficit 
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— 4. Corps géométriquement simples 


Dans cette configuration, on étudie l'écoulement autour de corps dont l'axe est de 
symétrie est dirigé dans la direction de l'écoulement (ou de la vitesse V). La 

— résultante R a donc ce méme axe de symétrie pour support, et se réduit alors à 
une traînée. 

Classification 

. Ondistingue deux catégories de corps : 

1) Les corps non profilés pour lequels la trainée de pression est propondérante. 
L'étude de l'écoulement autour de ces corps conduit à faire une nouvelle 
distinction : 

a) Les corps à arétes vives pour lesquels l'aréte fixe la position de la ligne de 
décollement, 

b) Les corps non profilés à courbure progressive (comme une sphère ou un 
cylindre) pour lesquels, selon la valeur du nombre de Reynolds, le point de 
décollement peut occuper plusieurs positions. 

2) Les corps profilés ou fuselés qui ont une trainée de frottement plus importante 
que la trainée de pression. 
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(2.9) 


E 3. Disque plan normal à la vitesse V 


» Pour 3000 s Re s 960 000 : 
Э Face amont : K, moyen = +0,83 


> Face aval : Кр moyen = -0,29 


— On en déduit alors le coefficient de traînée de pression : 
C, = 0,83 — (-0,29) = +1,12 = C, 
La surface de référence étant celle du disque plan. 


Comme le disque est mince, la trainée de frottement est nulle. 


3 Expérimentalement, on montre que la trainée globale est égale à 1,12. 
> Ceci montre que la traînée de frottement est effectivement nulle. 
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D. Cylindre circulaire à génératrices pe 


o erpendiculaires à V 
> 6.1. Ecoulement sans circulation de fluide parfait autour du cylindre 


. Larépartision des vitesses, autour du cylindre, est donnée par la relation : 
V = 2.V...|sin(0)| 


— L'origine de Ө est prise au point d'arrêt. 


> Le coefficient de pression locale est : 
`< K, =1- (VV)? 
Donc, la répartition de la pression devient ainsi : 


Kp = 1 - (VIV =} = 1 – A4.(sin(0)? 


P La traînée est nulle. 
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` > 6.2. Influence du nombre de Reynolds Re 
a) Etude expérimentale de la forme des lignes de courant : 
Les observations effectuées en soufflerie à l'aide de fumées, ainsi que la 
visualisation de l'écoulement au tunnel hydrodynamique, ont montré que 
l'aspect de l'écoulement varie en fonction du nombre de Reynolds Re. 
Re = V...D/v 
oü D est le diamètre du cylindre et v est la viscosité cinématique du fluide. 


— 50 < Ве < 2500 


— 2500 < 


амес 
2.105 < Re, < 5.105 
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b) Etude expérimentale du coefficient de traînée : 


plaque L а l'écoulement 
cylindre lisse А « ee yag 
cylindre rugueux Х« oo 


cylindre lisse А-5 


10 f 19 to* 0” 0“ wi 0° 
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p 6.3. Répartition de la pression autour du cylindre 
. Un petit orifice, percé dans la paroi du cylindre, établit la communication entre le 


fluide qui s'écoule et un manomètre. On relève la pression pour divers azimuts de 
> l'orifice. 


Ecoulement 
théorique 
trrotationnel 
Cette figure est obtenue en 
b> Ecoulement réel 


| portant les valeurs de K, sur la Re > Re, 
normale au cylindre en chaque | 7, réel 
— point. On a également porté la 
répartition théorique déduite de 
l'écoulement irrotationnel 


(parfait) : Кр = 1 - 4 sin*(0) 
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>» 7.3. Emploi de la sphère comme critère de turbulence 
La valeur du nombre de Reynolds critique Re, dépend de la turbulence du courant 


amont. Re, décroit quand la turbulence croít. La valeur de Re, pour une sphère 
`> peut donc servir à caractériser la turbulence dans une veine d'essais à air, comme 


la souflerie par exemple. 


On définit conventionnellement Re, par la valeur du nombre de Reynolds qui 
correspond à C, = 0,3. Re, varie de 150000 (soufflerie trés turbulente) à 380000 


` (trés bonne soufflerie). En vol, on trouve généralement : 
Re, - 385000 


Cette méthode n'est valable qu'en première approximation, car les résultats 


Q dépendent de la sphère et de l'appareillage utilisés lors des mesures. 
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>» 8. Cas d'une sphère 
8.1. Ecoulement aux faibles nombres de Reynolds. Loi de Stokes | 
Pour Re « 0.3, l'écoulement est laminaire. Stokes a calculé la trainée en appliquant 


l'équation différentielle de l'écoulement des fluides visqueux. Le résultat ainsi 


— obtenu est le suivant : 


— R, =3m7.u.V _ D 


oü u est la viscosité dynamique du fluide, et D le diamètre de la sphère. 


En passant au coefficient de trainée, on peut mettre ce résultat ainsi : 


. Cette formule présente un intérét pour l'étude de la chute des petites gouttelettes 


d'eau dans l'air (calcul de leurs diamètres à partir de la vitesse de chute), pour le 
calcul de la viscosité dynamique dans les viscosimetres à bille pour l'étude des 


Ñ  centrifugeuses et des séparateurs par gravité. 
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E 8.2. Ecoulement avec sillage et décollements 


Pour les plus grands nombres de Reynolds Re, on observe des phénomènes 
analogues à ceux constatés sur le cas du cylindre. Il existe deux régimes 


. d'écoulement, séparés par un nombre de Reynolds critique Re, voisin de З x 105. 


— angle 0 pouvant varier de 


* amont. C, vaut alors 0,176. 


l| - Le premier régime |: 
(Re<Re,) correspond аи KS 
— décollement d'une couche N 
_ limite laminaire pour 0=85°. SS 
On a donc : С,=0,44 AM 
à à La de 
~- Le second régime Stokes 


(Re>Re,) correspond au b. Wo _ 
BE 


= 24/Re 


décollement d’une couche 
limite turbulente pour un 


110 à 140° selon la 
turbulence du courant | -2 


Deux disques l'un derrière l'autre : L = distance des 
disques ; coefficient C, global pour l'ensemble des disques 


Plaque rectangulaire de longueur L, normale à la direction 
de l'écoulement 


Cylindre circulaire de longueur L et d'axe nomal à la 
direction de l'écoulement 


Cylindre circulaire de longueur L et d'axe paralléle à la 
direction de l'écoulement 


Hémisphère creux : 
- Concavité vers l'amont de l'écoulement 
- Concavité vers l'aval de l'écoulement 


Ballon dirigeable 
Profil d'aile 


10? à 105 
10? à 105 
103 à 105 
10? à 105 
> 5 105 
> 5 105 
> 5 105 
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